


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 


DSpace Repository 


Theses and Dissertations 


1989 


1. Thesis and Dissertation Collection, all items 


Why stochastic modeling is essential in 
analyzing interdicted transportation network performance 


Chol, Byung Kwon 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/27335 


Copyright is reserved by the copyright owner. 


Downloaded from NPS Archive: Calhoun 


atthe DUDLEY 


WN] | ciseaRy 


http://www.nps.edu/library 





Calhoun is the Naval Postgraduate School's public access digital repository for 

research materials and institutional publications created by the NPS community. 

Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholarly author. 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 












































































































































































































































































































































































































































































































































































































































































































































4 ‘ . ' * . 5 - all '* € »* _~ ny a 
‘ ’ a4 i Ce Jel iy ! x Vey! eae « BB! at ‘ OE FG ELS TORR Ag CIPO A 0FMs OF Ne eget (AD, Aas EtA Mle Me atAl Alar Rsk: Cig B7e O. a7dqt 1g Beesley C48 ON 
s . ar 14 ’ oe | te? my 3 ‘ ‘ x . OA. FO aed be EN 1 MEM Ble BW dh bet NO A we ubanoginns % aseondyeit A Mad Atte Beatie) bi achaneet: acat-h anauenty 
‘ : a) re ' 3 wee AY Bt I *, ae Leeds ERE MASS eee ERR AEG. lle hy Ay ed Mechs 11g Ace las 0° 8589 6 11682 IAIAT OI Wi ASA, 14 Me pensar eras 
: ’ a u ; we eer SO ad CREATE HoP AW AMAA rt BGO R SM Aid be & ade Kn peter Mem WDE As as iSksanss mateo ee 
' ‘ 1 4 1 ® 3 : ae ‘a Te YG 8a pase: eset eee redo Mirtle Layo jor A BALE AI MeN TAp Me BOM ALIA, TAAL MASP LOO Bk RPALATNE os WAL MORSE SEAS Oe HB Oates 
; & es “ 4 Up 4 Ppt ot PAZ UENO ge 6 Uk Rar) berm PA edSut@ date. gar, ages PAM ae CO Doe ake Ce Vy Teh AF et 2 OMAR RIAA Miao AhraWegiere 8 © aaeest “S10 fas Oe Sab eat hes SOR ve 
Ly . . Cone mat, 1 ] t sa aisas | te at es ' Ie 6A 4 ot Nhe ne Maite nied Mae Cbd fete Se Fk MAS RATAN hg 4g OO ce a ated, te Pelgu gy a Rede BODO TL AG ALAM, fos An Rs Be OS a hs Rete rs CometgD 
mot ’ ‘ e 1 ; La a o1 %)& ‘ oasrging @ cept 4 RA La Wen) Bagh tree BM Epa kts GANAH. by bs B90 stihl COA Woh TeD 6a Ae 10? ae he DIR LA 8 4RS ae OE AO SW Ree, 
1 . a ts ‘ 1 i eriag raat a‘ BUA WeYtsts & ty ee et Lied LET [emg | = peek. ee Ae ae a HAP AL BOA he fenton Ane | DL te! a tte A541 be PHO ALO M6 AQel Fe4 7 80: ERR Te 
, - 1 1 » seta ! @ avo hrAt Rite MAE '& Lotte 041 bo COU SMe Tbe, BaNe le UD OURO rag Me As ind CRA AER REIL Bar, PQ BRAD ph LAL AL UTNE aes Eh OL Meduuw It Bs quer a. W Seqy Ga been 4s agen 
. oa A rs + irae UM iced gg 8 BOT oN ALA BEA Leben Finda. TRA AW ORAL, Ola shit ga Ae rag ht 0/8, Bee ghteM, 0.8 adtiat 48. Ta ONS AAO Bry QO. snd eh, OF Brads dy Me W Aer te ye Of iemenee 
] 1 1 L#, Proll’ take 2 Fo oe OR Oot Olid O SFG we Bate WINE Mra (ahaa aha O0e VAS Ohateo cay Asa. e. he a6 asa nd tS 0 RL a be AOC Og ty 081 he TERe Oe 
3 : ‘ 1 6 re wees on 11.0 Aly e Ad Oust oa! ott tse en ee eee) paca tiga W ptt! { CIR a Temes UC RAAT Os Oe AOE AAAs ee AAA, EN Dy Meg BINT Liye Oe BARTH BOT R, Pe 
4 A : ' ‘ Co yo COT i TP ee Weta bic a as ee Ll rg Pr Tad aoticnisad vacrvonty BA LAs Ree 8 Fh Oe O50 As este & AV be We Nee ee BeReCS 
5 4 ‘ ' ! ‘ i ee gaa A F ee os *. toa U9 tay % Leys eR Pe 8 WN arty 8? Riteay Saad. PL te ee ho ie | Fie td Fa en Sah had, FO DAS Ty tts MY Rd ROR Rm 
bd ' * a‘ Auk y SUT 54 Aogeut “ Noga Gale ats f . SMeetIOt A Griegates 1/10 Ava im Rae et > Qe eB dh Pre er 
i a hy ' ' Sate ce Cas Me a ® Wh ye Fl she ty ae Sie Ol ey toh t 8 aa, ibe Ms gs UF Rates RTA MS Mary 18 FT Ol of Re Rete ee RM FRC SMER, 
i 1 » ieee t ts + bh ide i 6 ‘ ] wave We enlace AS SAG Zoeken ays "ARA@ {ar aon et t eae Rat! OF adits ha OH ee OVE AMA ese Wan PERE 
7 e a we ae 4 i] dyn! > Usa rae ’ Aa Ort eR Mohs Cine PROT MI ete ee Pot eh cee Le trom 
pe " ot phe ' ‘ a ' 1 A ee | e* a0? © av abated Fou Uy ='ah of palar Vasaae area Or. Maryn & 11 Oe GURY OF, dy 1 Dp eg de Oe 1 BAN OED 8: 
uo: 4 Lane ' Ce an 19a? , BOG hh 86s Mot tet Aide GAPda U9 WA ® q XA pt ha! md 7g V9 yD By Ohl be UO Ot er Ad Dy MIG 18 sear MANs By Ri Cy RD My debe da toh & 
; 1 a * as a) re oe J' per raerah ee ee Vat erga we Wa, Lue, Om Me ad WOE Le hal OE ANA Bee my oy, 2494 1 Been Ey MAM Aiele Bg! 
' « ° o 4 a ri a : 1 0 OPA ae we 8 Bit stesG st qaQf 8 Nore Ju ta De CMA Lae A GROPP IVA see Gale Aye er Pe ee dei tgh my 
: Fs . 7 , ' Sor hd obd ey AUP Ts Meheh. Whe Tanetr@ tte Ae @iar WdeRee Agha cs Moy ay fils heute oe WANES Fete Oy ey Uta an Os Onto 
A e wasn es r : r t« yh e'g' 9a P&L . SAeesit fay tt aye, oe aru Reaf@ wp Hie? a 
' me P| \ 4 4 4 ys a aL 4 ; al ' 7 , ar " a HA “Hh rch apn A aa Tg rato wy OO Vator ve vs 
. - : ; 
. ‘ -, = a eR gaa gotta tin feng 
v 4 ' 4 “a 1 Os ve uw S37 afte 2 ot Use LAs Oa 1 Oh ete he Oe OF 
. - 1 ' 5 aia + z ‘ re Rg rhs nL MPO PRD Ae 1AM Or RS AN 
; ri are 1 % et ° ot; ete aie & a fed vol, 000.9 . Vg Meee po igen s Ot Be Anes 
; a? ‘ gy tae feat gag ae Uh peek s } she 1 Winter DLR, No Su LAD AA Powe Fa gtah gO BUOCae ASRONAN SBT! 
* ' en oe es ee babel 01 gt Sra Pad eis CARS Oe -9 re 
' 4 ’ ‘ eae 4 SHORFARAL IY Rasa, A eS Oe BR Reds 
e ¢ Sean on? cae} i Oe a whim one & Maeseanas 
. d it we ro BA 9 Rede ie, 14 %ef mah 
. ‘ : 1 it ya a % s ut wt - : i. Ve 14505 my mM. cete 8 IMs Ul Rar Me eI oP Lies 9M gteg a 
“ Ce atk iaeer ai 9 ws 1 . sate Ore te M eT Tt y 4 RINE SE OT das Fl Range Aaeargaaeae cane hs Oe 99:90 th tne 
+ teas *o4 “ mer ° i eo pM AT BD Arle tate ny & OE d hy oC FA ie deem, ‘we Ge Meret sate BPM Reee.% Oamine 
- ‘ ae ‘ 1 = 1 ' es .° 8 rT] 14 1 90o 669 1. we ertsbhg Fe 7 ‘tetra UHRA @ favahewany fps onal ALP) fe OD aR an ay A) Be ¥ 1k Ce, 
F Shan 5 A ‘ ar par Pe * aN ' @ he nega es eu,%te § tyr we LUFS Dan IAD RP Ady Vata, Ia hg ki 00g AS Py Sy ee ag 
+ * 1 e a0 ¢ “bs U 1 LJ s a yh WAG MVE ses BUA EAs cre Aer As824ons Packailatiee aa 
An ’ ‘As ates : ; rs F aunt Par rT i re « «Va ats , Sete de on a 9% Siem te tsar ta8.8, 
= A 5 Py i an a 1 ' ! cl] 4 y 1 wh? sto Mal * s04atnay ror Maegay ange 
, . os ee q aig Commierait us x ye 8 ES AA 8 Re tdNIAGN & ye hiiotadk tT abpae 
5 x . ' i] 4 q A AM 5 mast deme ls 
CRT ain en © ‘ e 1 ot ee i 28084 A,het 8 yay ‘ tes 6 by By he i? 
’ won > Fitted o a MSE a brerps te Rengermateg 
r Fi Py 1 1 13¢o a'se te Orialn | ham Staree 
; ' > ' Da oa tM gt - PEA OLMM Ardy Det 
5 otieet 8 " v8 1¢ z 325 we Mee a: saavarn seAyhthe® 2 Ps May Sewer e wagi es! 
' 1 = 1 : fet c; ee ch 1glAsn eo adit 9 $OP a» yt dels whaeakie Hat 8M Res Vie te LENG ARS fy ey EAR BIH ee 7 
7 F ; a er , : fro uot > at peat tk Fe 3 ems 049 Geept eo Att a! PaNeD ede haw oh ete g Mtg’ 2 RP HLA omy, RIM o NEM. egtty AD Meat Rey eRe anaes sem 
; s] aa} ‘ aiid 1 Eat ni - wes ry 4 fear ign 3h » 4 PUVA eM, WATE MUL OMIA Oe eae TA Re and Eras: oe DATE em Ogre Py dy inh ADTs Heme 190 he UND 
. cans or, | ' ' i as i ‘sic ‘ } SEAN @ PW hi Ayerg sree pes eM eet he ays He tees ie wat; MEAS 9 hr Ng he TEs AL bs QI Rae ye sO) Fe 
' 5 ' a er} ? ay es | Agata og ‘ AA Ate & O18 HdtE Se Spe We OFS, GEgt tdubo 4 me 400 A De wpe aa baal 
: ‘eats ; F = ' or iia ee * ey US PRP 0 Mies Ghee MIE ee: Om be ccs mrasendd “Moral : ee es ache 
° ’ > a1 Mls 3S Se eo Wale yrqcemae Leen "my Marer Cate tte beset Oe el alge 
3 id F ce r a. 1 ° a‘afete 1 i wari © aga WSagryet rity Reiimemara,! pee te gare rag l, 
‘ A Lu , : , ‘ se » * we pig 8 metre! Rens tone n ee aM Rote bh suey ADOGRA Tents Ape Bereta 
: : . : ; ‘ ate ‘ aren | SRaW whee c t.092 we Shen F pease sO 4h Perr te 1FDF BRAG Ae & FOS, LPI 
re ® ' , ' Ue 1 J Lee) © TA Leh 8 Og bet og FAS Be WG FAL Ep TOhs CTs VS Ag OU Lore harden e, Me Ronde fing ug: Mitr: 
: A 4 1 cart : « « tar o t 1 ' a Pat ay ber a Pavey mo ee Aa! te Agey Ral te Rem te Dye ae Rett alos 's (eNga, oa 
F _ : wove 1 \ Tye 7A i eh ee a oe Mm he ne a | Fie aes (* ted nena, 
: » : ava > A "@ 6 * , 9 8 pe te Di fasr “a * & herrhy Te 
' 4 20 aout ' " fy 1 ot! ' : PIAS 0 hime G O34 ; wh 
? » t wr 1 ee ' 1 s4 os #6 Q iguder * 
1 4 oa 4 ' & i y 1h Pept ees Ande th "gaa Uni a 
i 1‘ 1 atte fel Bee geetd OWs hee rte ¥ 
A ® + 7: eae, 8 Ph? Beads «pease ® 2 i adene, OP 
ry - . ' ® LU Py fy a Mise pare 2 9 hips & 4 daly pe utd aetna our neeabeng 
i 5 A . U 4 & z fe % ZIT Te dees hy vs ea nas gin gs ae deere SA) In \ et i hee 
' ' “4 , ' pate fey 4 a ea . * © dF eh “tek phos ofa go htteoe Dt ted ed Ag : Oe be ab Fo 0. An WIM 
a {ies ai. it he ee ae es Cady Da ems Eaes met 
® oe if t 14 Get Uda? ef arn rag 
5 » ’ ' wh: 1 ' 2s 1: oo Deb gge is lg: ati eo tidgloss AP 
' ; as Pat) a a 4 ar ee eee ee Oe Sc et ee | Pg henga erat - 0 sagt Of 
a yy 14 a a) 5 1 YA DO" eae e kre te US PSOE TR Ae Mn tote “Fgany, nntyae ate, Ra RgPk Deeese re, 
1 : ' ‘ ss! ae * Bee he iat TaN igetg, 1b omer lhe OFeutipe  Veasyreo sites eye h UMredass 
Fi = . Sy a Hy i one i] ray ent okt thie Teun pega, AVY fee cpt en Hat Wet AgRi ey 
1 s Ms %, 8 Beem er eeanges “Ady 6 Neg tg 4 
' ' y % ‘ . & . % Sae he % 
’ - a ‘ a] PI ee. a ' ie | a@e ° he ¥ eng tow eee ’ Aparsd 
; ry Ee one . 15 J ater 0 a8 Boag ryemqe Sane ‘odes z 
° sere i as @ : TM Yu a ON Spates . ‘ "nef §o% go <1teee 
r ' his ' Ww ‘ e 9 toe we 4 * of 4 Fa 
a e 2 ha 4 : ® . er . 1 Le swt Ve @ at . Ry® Be Vgh . Sage 
1. iu a 1 ‘ 1 ee sus 4 ' L e: Tel WA wr Anions 4 
bie ’ 53 =o . ae | ae tpt 1 i pees 1 
1 ' : aa . : 4 ay » ap ye ,! tite 0 tee em oabetie 
ws ae nies ‘ s eo aatite . oe we free Va oy LL ed eee eek it CE. ve Pwr err Tity 
= . : e 1 ‘ ‘ P , mn at ae 71 gu 8 = Sy, he HELLA Osh Ar ERAS ten oh KUO S eer hwe how ol gra te be oars osgn farne ahesos ona 
4 1 ’ ye * oy 1 a sg. 8 8 , 08 er) n at %es pod n tana’ Sa v8 ag een ass a wm! agin ta St en ee ee 
2: . . a 8 8 sgt O49 my fee & @ 1th 88 yore ay eee rr «a ee f . 
a . a » 1 ak mo ’ % s Be 68 Medsad ty - 1a4q tge Ame? tare mimi oo Me Pedy ® ge Qe @ BETZ U rales. Nom” phe 
. ’ whe © ‘ * ‘ . 2 uie a Wi oR Greet tae t af . Apestie BF oy", fe mose oto Fos amazaate 
o Ue ° . y\s »ee sat ng pete ay oe Ae pp UP ote MULES Rinst ete te Boe Vem sy ce, aye Ae, wr ar bv yo By arene 8 
a % vad) . : tes ee ae . ae alt =n 8 . + dap mphie hi twel ache ha cy: whe . re eTeont nit Seqen 25 40 /4~ 
e Py ot 1 . rr iets en ee ~~ Cy PMA Boson TES, BOrhy tpefeye ab, Sikes . fee Ss OFar yi Nut Sp Rese pRmemene ter 
1 * . ¥ 4 2 2 44a easy “$2 Rose i ty Gee 4 1, Grey t,t dere Ato. s 18h im berbpia® gal" Pays fey Bee a 21 4 mie me Oe ees Me 
= . 4 ' . ' Rhrrey er Po Cre | s+ my Oa ate Ae! ae Se 8 '~we ws eh B ts Ta. e4 a teh ety ee 
. * . " te Vis ode 1 p the te 0 ST f mae) tee tan en mone ay eete & fs ae vedas 
' 5 . , 3 4 a) es hue * a tunen tre! Che fese an, “hte Tease 20, Deseede 5 
' ay ' 2 2 Sheol » tnt > ee oe Ast shes el bf sind 
a é * ‘ ry . ° m™ > > ee @ age eee fas sifce omens a 
; F fi i ASPB EVAN, Doe me % * : aad clyeanh 
1 . 1 2% * % » ’ a | > re Pe bf op ha Sp 3 or alr s Wiese yal be tas 
. ty Z ' 4 ] Y omrdaeds aly be» FA teahy 
1 : ’ ‘ . . o © §Ms % er yt Yess » fad age 
- 1 Sizts aofe “9 
a? e : ; § ' 5 ae ed 
» » “8 ry . ae me teh Aisin te 
»! a ae : : fod ' ae oa iad Bas 
Ly i : q a ade 
voor See ._  & igs ON 
o s f le ‘ 
. ‘ 4 ° 
Fi ' . fet 1 : 
. : : ¥ \ 1 4 * det 
LU ue “ , Ae * tag & Sifoinge 
D -° . f tt te tM, Wek noes 
ts a . a} 8 lw » lhe , . wai foftne a a-ts seag 
t : +3 5 5 ye 77 Py *. hse & te FA tr, Ped Thales of e& 
® P ir =e ' As eo! t dae Oh vis toe he fey wtp nth 0s? "Sph 
. t ‘ qt te ie Pe 7h ") » 
1 1 ” ° L) ® . . ae “7 ee Ls as taat = ie Wu om! tarts b> fies 5 
‘ 1% Mie i a 7 ae Pat 1S WF fap ot "Ks BRS UL ag esa S pret Saeeh as " Are eet 
Py we Lé ; Sea sietis “aloe , | pe Auteerte Sl 4 
, 2 an efase SERS 6 ek ined aa 
i “ ® * a erate Lite sheist ee nett: a ty 
. 1 ° 5 ’ Weare neers ties Gs dhe tt wag, a 
: ° ’ ry # wt, i : sere Beate 
7 ® ‘ » 7 area Pura. 2° fe) nga Fae FN ody SRE 
. ' \ ra) ‘ us a? . > ' % a” @ at Sebati t. ma thy ime ; v eeGhitahy f747 
= 4 r] yy 6 ft i ¢ P “poh : r +> i Rh lit ee At? thx? 2% 
ae a ’ * Janae Us i a ee a, B PE soy 0th 
Po ie hed f 1 oT ) ' tert 
, 1 . es . Come ad e e 4 ed area ' ou 
4 i £ 1 & ® ° a2 i ’ 
Des “4 1 e ° ‘are 7 (Cres e we “make 
e ' pee ' " 1 . 0 ' . ep ‘ or 
' é e . Ld + “~! 1 of ry i 
' 1 ty ‘ 3 5 
te P . ., a8 4 #0 dl ’ 
mane . a 5 +0 : uJ > ik , gf 1 
. 5 ® 1 ; Ng ery # ott gt Mh er ow a 
. ¢ o- oo 4 Ce FT aes | iA 
cA . wile . ' ' . ee oa — ste ae? oh, | ty ere anter Wyseerereses FRs 
. ot 1 oe. v6. ca: fis wwe 5h ow > ae 14 Ser emt a Fk Site wees ee AL 
*, Py ? : F = | Le eer Ci hae Auge 0 € sq 7} 1 LPTs Ss We os tah i) heed Cd 
® . 4 ' mY i Pd a | t k s ier paps Se 
5 e ' rs 4 1 ' . 1 0 Me ff Berets % "72% - Sgt eet A 
s ' 2 oe, e © etd, £ oye ° THe Zc Fert Ut. 
se 4 ‘ Ly ‘ 4 ts i] asl fa 
4 @ . ¢ ee! ss 
5 ‘ , ea o Gay 1 4 ' te) a a ee A aenk, 
Py . . e v ' = 1 "4 t' 4 eecea 
‘ 7 pi @a 0 et e a r] « oft ‘ Chae! Car iy . r) 
. . 4 Oe ll ¥ } ay ~ rr 
: . ea ; ‘ * " tee .? 
= ° 1 x alt a Sie : 
e e ‘ . , . L td Zs . > 4 f ? of" * 
2 yee at ees Pa ° Q = 
A - . 4 * eo. ' ‘a " rf a5 Sincane i¢ + : 
é 1 « f , y je Le yt note, at 
eu ‘1 , ry 2 s . . le sae we ’ 4 « ' noe 
1 ' 1 , ens * J of yore rd ot ed eh hd lad 
‘ eget ° ee s ° ey 7. Oat es a! ot" tte # ot IRS Ra Peete tee te 
. ' a , . 25 x Lal SS © mee eer ee 
‘ 1 ° ers ot is @ oe pre eMmre ‘* 3° * oo Ree e wh Ne ewarn ents 
s ' . sare ' FB Ons dg ae Fee ee ewe ye Oe eT 
1 3 ' ts a 1 a 40 #8 «4 re Oe * te g* fee 
1 z ye Fy 2 . we PY La 4 i 
F , sof . “° “4 wou asi et a a eineate CAC ANE rte ots, 
' . 1 . ' “ vu w« r . cae) wae As TL tos Ya tool Lt a dealt Ll SI aad ol a ye : ” ana ove eum: ° 
. os ay fas ee Fe sue t ot Cha) ae) sae glade BS (en ees ee. YY Ale sl aesere : ee SP Oe Seles oti egcewe oS SgkOn Mo PeRD F leer crarshd eotede thse ee pit aig es uaplps ot ge wae choceushgeceete 
4 . .* ¥ ta of Pee oe Fore ty write Gg 88 6 Otte Fie og Vere fe RAITT FE oh ge Net Oe 4 tates ae el cot eet tg 8M 08 ob Vey re" gets vee tor te pore Ae 
eo. ' a . s ot ¢# » ' sd ee) ew 1st ae Ite sie SoBe gine myers © Ogee y Eth ® 1% eons tyes anf erBts ve ehh ete © precy sae st : lis sammie ge 
’ . ‘ ? 7 of #8 760 wre seit ¢ Sr * YAe soho ee Oh sgt © ate: UF start g seer sesh, (oh etm aecyetas PALM 2 fn-s0"o 08 
' ac - wane ee ' ' chr a Wierd oe é set, Pere rwtrrcere A eae ret ef rtases: Ped ERs agile a AS go gy gates mer Sprerrs Ke 
- oe. ou ¥ e ’ ' » we et 5 6 ete Ia are stat of2529 tt SP yponsen pe wa ges A gh nn lee A gt gee pe rosa ee 
: , vee as . eal ' - 8 Las By LORE AS Soe eae P.O F 2 U8 wo Bia, sqthgbernt tight Gem enThen pie es ot wares ss anene wee 
7? z . ? é 900 08) a : Me ,e « ore WP seston fue pert Lat ee 1 de nee Pe of 8 gn! anes 4 
| Pi ' 1 .? . ° ‘ i oe stv @1fe oA He - FRE Ra eer gaps? LE IOI OF gs gH Fie, Sag eee te eer eer* #5 
1 te - 4 ‘ 1 1 60 098% 20 PP aes "Sen Reve Nets t <2 arte Lipsi Nees pass ¥ ph bp feat es a 
, : = e e ’ ® . ' . « > Fe ee Pe dedi Me? Peleg Cie op ab agh veep ot anger Set ens eye RE Oe 
. : ' . - ' ts il F Set ee oeieate Lad aT es Pegres 1g jog e Sph ett yw Ua A = swe e.gurenen “Se re fe 
e =n o ¢ - @ mel ' ' oe? « ia orm yall aU pec OI ct oe “pens Cap ee be | ” We urcpse mp ovpwnry 9 Fiz Pr ee 
* . ; ' 1 ¢ wit 2 Lad beers GET 28 ye meee! ergs ote a et ee ree Oe 
1 . . ‘ 18 ye , Je es 7 i ’ ' on ea 5oeer fe 2-7 # Me saerne 6 © PYRE ey trs gre Fe ee Es wy 8 
tate yO ; i it 1 cin ' fu “gt FRer™s re Nad SAO tespcgeetyte- ce Sit gs girl? tbe lapels ee WF * fe 
onsen 1 . 2 sr 6 ee oy eZ A ene of @ : - 0 PR otek ey TPO ae se” Bn ge TNE Teer! eRe 
: Bs a ue i : 5 ° oh vs ' ’ é ' y BOOS E Wiergt Ss Of eG mt et ot mt yma Fel gine we! 
' se | ° ee tes mu ' Put © «2 ge oe Fiesty bs , FR feo 29, Ewe hon Yer Ore! £0) ect a? = ee 
s i A s . chs a tes te or r® «Per gtig, ee gcisaealg Tat oad te as lp sat eat an mb an FPm er Grey 2G) Cyt He 
' oY y at , oon % 8 ‘ a a oerey ° ween, Fn erat, we lao Oy ae a Oy WMY Reet ok 
ae F a « 1 le . ies F 1 sf Pi . o@ Porn pew ere ye - ere ead « 
’ oor’ . 1 1 Pat bet) a. e ees 4 1 te “ . eae . a awe ie eos PPR sratare © 2 Fw ' ow Wm afta . 
‘ Fs A 1 ; v . 1 a o br ow ot ttet hte 5 te S0 tye ctlee orytiss stan 0 ytonwe 
. al 1 a . cy . #1 Pe ae | ° r) ot a2 @2 08 ar et ge PS On h Oar (UTE! tne Mots 
‘i t ¥ 5 A iu 1 ais r . . . ' . os bo Pn eT a et! peel + ol emgs fg f- ees yatnaesiis ‘ ; gatas went pen 
. . . e . 1 ' oe ee or | ' g wa € oe ft 8 oe ke . re PL a ee. a eee 2 2f beg tleats Rsig= Rl FOF eee, 
1 . 1 ee 1 y 8, 2 @ee#s Lf vA .eite 90a Me He uo TW SYOly seqsres ae pete pron 
. 0 ‘ 14 4 . ‘ tne ‘am 8 1 1 Ae 1 15% ® ony . Cd o “te ‘whatqa ¢ greet 7 et star de, L0aee He, em eey 
‘ 1 ‘ . 6 ' ’ v ' ¥ * yi at onee Fern 2 We 8 Oe he ig Pr td Aw, sat oe A OPE ES > 
: ‘ Fr @ 4A wae ars ' at Por) at weet On a a ee i Pyrie 0 ohe hp ete Sete, wrens ee PPI ese Pelee se © 08 Qrter ye ter. oF MINS TOME nET ant 2! OP ene eer Ce of oh ye 
a 4 a, ¥ ' o “ve 8 82 1 fz ‘acs C# ‘ ' Sed re Fe dar aeas gram ag EQ OTR' 4 a esc PAY oe 28 ote uk LAD AN. AO PL AT 1p ARE PN IT Oe Oe ape ee we lee ph wk BTR eed 
r < - rear 7 ‘ ? ' wee eenwm ene Be Pay, an 1 owe ' eogter vray: Rte Sah eee errs eet vemer eyed 
' ¢ it ‘ e e n Char ._ Pa nd 121 6ye oe gee Sar gm iter 18 08 wht eae et pies OP gE HY e yt) 
os eret iu oe 1 . ve 1 ry as es 90 1 ' . Wott re Wl arse, 4 ep Te 6 ¥ & ~ = sree cQwhuss 1 (fe eyes 
‘ rae Lng mie ' t Fi y wae I 2 os Hot, % Apis & oe mPpeVEes wt rete 
- - , . g« e ’ : ots ri ‘ ay @ Ooty . ee Ws FO > FE Glee ¥? ole 4e 
1 ' y te 1 P _ ; ay ae ae £ af 
. i s “ ‘ a e ae e “uo 
Secl ge wilt , , . seus re as ; a ee 5 dit he ay. oo) « tt ot te othee ss ae A AL NEL 
4 me <n ga op ae $ oy a+ fF 18 46 feet ute os Fa AU ga I IM pf ipeedine “% 
Fi ; 1 . * Te ee ey " 1 ¢ eta ewe At we CAD erp = reyer es 
' aoe . a = as ra) 16 a . aos Y iewe Pi . dea LJ $tebhce eae wey 5 peer ete ghz" 
er ies fe le ere ees 4 ths WS ar td og phieleerne Wt 
» raat £ wo: F aye ot on @ wn bee fe un tee Mone £ (he of etewvrmrTel sce 
1 . r ' e 4 ue ao. y Le r- Yerherete Fie yeas 
' naar * a te » 2 4 1 ert eo 2 * 8 , : i a ee 0 3Per 1 we yGee Ite ge oeds « 
ee : a; ‘ - Fe eo eat HP. ps | See ges Pier hy 
oe a ' i] . le aeptee het te | ecreMe sn O88 
® F ' ra berg ' 1 ‘ ' -* ef at 
= aed ; 5 1 nN os! s 1g 4 ' . - bee . 
= ° ’ he ‘e g 
. « ' ! ae 7 - 
s ef ‘ Wee » at RTeeat at bean 
. rod 1 bd aT Cum m7 e Le t oes ' ie pg He LM eG 
’ ee » 8 Yr 1 , o” ? ae 
- 1 ' ay a it) , i a ey ¢ o Ate se nO 8 = . 
c eae he a ere peat ye aes, ar OR EE VTA FD A FOE IS A NTT PEGE He 30 Oe WOTTON PS TOK aes Dales sim e'e Viglern eh ener, 
s vet wes s 1 * a a “ene os wy zi ? * . i J EY Lhd - ©; 
: . * : e me : = a ee ae : ' . es 2 re ' ane whi ret 1a bela FU Ne call sentry siete een piglet pri potas cme Ruewurre witene govern 1 pacaepanacnedabe eel apap le csc at ie 
ae ee ; Pra eet a Ta aie = . if re te me Pe © eases. TE ie OS oe Le ORM Breese poeta per eee ah o 10h of gay rate MYTH oratye Pra ra ye pig hte a ae 9 Ee Vy, . 
P ’ Sialeaiin nae ot ee oko hl ’ ' - pt aE 5 i Akar oh OA ph Ab ceeh a ila Papakea el toate nal t 
r ane Po ari rn v4 nie : : e Pfu rep} fs ore Re ee ae ae el a ee 
4 ve C eo. of via c de £ 
’ 1 e ° “1 oe. ye ‘ . ew. ‘ot ? 4% 1 i 4 Sled 
é . . ° of . - vane 
¢ ott 1 1 ’ 1 « 1 Ve * Ant 
1s i . ° ; " ‘* Dit e1 vg tee i »er . sen nt) o [eine 9 195 He OB) Few ig AY Ue dl 
1 5 re me Fae ee 4 i 13 fest ws 5) on Me Le ite TA Serer A ANS HO ny oyeLEn stderr e acy: 
; a 1 “ws - offs av s ae & ies «0 wag vs tts. eye wires ¥ yore wee « 
1? . nT Ua Tey ie . . = Sent * ' a ee ey pak au # Aye A pe arn 
tele ena et 1 ae ® ’ ' 
r é ‘. . a | ' on 
ry cio ot 4 an ary t . wen 
d 1 ott by bio “8 c s ° Pulp vais 
a ¢ os Smee we 1.008 6 ¢i y i ell 
1 . a wi 7 ph F 
1 4 we os si 5 2 BSyore avhe ose OMe 8 BUCS FR i 
' ea pean ames) okays ono anna 4 ie Me res on eAcwe het es at meatal sunli be ante be ee 
2 8 oer yu 194% aa r a* yee et @te. tae ais aee ty 240 Pe wise aes 
“es . . 7 . ' ' hes é ¥ @a 78 a ie a a i 
’ ’ bie (he nie 4 ‘ tm pe ity ete pe ee egy age 
‘ . 1 . . ° % - 1s 9 fe tie Fil he) pee be Ll mh ¥ 
of 4 ; . is telus dee U a <¥ :o1 ce an! A Da Mote eh ello! 
F .« = F RS fo. et a NU iE ioe BR Dt dk) 
i] i] of eo? . ° F, ore 2 Sift a) ¥ Sa 0 PTNNs FF O° gute 
, -1 ' . ' . 18 ‘ rae Pp A babs] 2 aD 
1 . er) od * » oo ae TEE ae Tg es 
* e ° ‘ te ’ 46 i ert 2 ~: ~ e: ~ 
i . a Ye [a é ¥ 1 ry a re eo 4 
' ‘ 2 ; ‘ . 7 st Le —s de te x o ert We 2 
' ‘ez ce os | b «4 . rt eo ¢€ ous = waa ve ve 
= a 5 ° ' 9 mt é SOT wt ne. er 19 tL Ar ed 
‘ » £ e 1¢ OC yank ey rs beens ee sa¥. OS West 
, 1 rot « 6e vets ss mvs 
5 , ; F : ae Pe er te lete ‘ Sr ie We te ted cS 
a i 1 Sen ‘ é > epmae ahs, erm, t¢ 4 ge iJ 
‘ ° ma kee scuetase ry $4 side Geog she te FT MF 1is09 : we ofa Stee ras ae Id os 




















NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





WHY STOCHASTIC MODELING IS ESSENTIAL 
IN ANALYZING INTERDICTED TRANSPORTATION 
NETWORK PERFORMANCE 


by 


Choi, Byung Kwon 
«*’ June 1989 


Thesis Advisor : Michael P. Bailey 





Approved for public release; distribution is unlimited 





UNCLASSIFIED 
SECURITY CLASSIFICATION OF THIS PAGE 


Form roved 
REPORT DOCUMENTATION PAGE Gp ielals SAGAR 


1a REPORT SECURITY CLASSIFICATION 1b RESTRICTIVE MARKINGS 
® Unclassified 


f 2a SECURITY CLASSIFICATION AUTHORITY 3 DISTRIBUTION / AVAILABILITY OF REPORT 
Approved for public release; 
distribution is unlimited 






# 2b DECLASSIFICATION /DOWNGRADING SCHEDULE 








{| 4 PERFORMING ORGANIZATION REPORT NUMBER(S) 5 MONITORING ORGANIZATION REPORT NUMBER(S) 


6a NAME OF PERFORMING ORGANIZATION 


6b OFFICE SYMBOL 
(if applicable) 


Naval Postgraduate School Code 54 Naval Postgraduate School 
6c. ADDRESS (City, State, and ZIP Code) Jo ADDRESS (City, State, and ZIP Code) 


7a NAME OF MONITORING ORGANIZATION 





Monterey, CA 93943-5000 Monterey, CA 93943-5000 


| 8a NAME OF FUNDING /SPONSORING SbeOFFICE-SYMBOL 9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 
ORGANIZATION (if applicable) 





| 8c. ADDRESS (City, State, and ZiP Code) 10 SOURCE OF FUNDING NUMBERS 


PROGRAM PROJECT 
ELEMENT NO NO 









TASK WORK UNIT 
NO ACCESSION NO 









11. TITLE (include Security Classification) 
WHY STOCHASTIC MODELING IS ESSENTIAL IN ANALYZING INTERDICTED TRANSPORTATION 
NETWORK PERFORMANCE 

12. PERSONAL AUTHOR(S) 
Choi, Byung Kwon 


13a TYPE OF REPORT 13b TIME COVERED 14 DATE OF REPORT (Year, Month, Day) [1S PAGE COUNT 
aetewm Ss Thesis i re | June 1989 70 


16 SUPPLEMENTARY NOTATION 
The views expressed in this thesis are those of the author and do not reflect the 
official policy or position of the Department of Defense or the U.S. Government. 


7 COSATI CODES 18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 
FIELD SUB-GROUP 


19 ABSTRACT (Continue on reverse if necessary and identify by block number) 








Road Network, Interdiction, Transportation 





In this thesis, we use simulation to study the probabilistic dynamics of a road 
transportation network when the network is subjected to interdiction by aircraft. 
We consider several road interdiction schemes. The simulation results are compared 
to those obtained by using deterministic expected value optimization techniques. 
This latter approach has been used by other researchers in studies of interdictions 
of networks. We demonstrate, through the use of two realistic military examples, 
that the deterministic results poorly predict the performance of the stochastic 
systems and that the bias incurred by using deterministic methods is significant. 


Therefore, the stochastic model should be used in the real situation. 





20 DISTRIBUTION ‘AVAILABILITY OF ABSTRACT 21 ABSTRACT SECURITY CLASSIFICATION 
EX} UNCLASSIFIED/UNLIMITED [[] SAME AS RPT (J DTIC USERS Unclassified 
22a NAME OF RESPONSIBLE INDIVIDUAL 22b TELEPHONE (include Area Code) 
Michael P. Bailey (408) 646-2085 
DD Form 1473, JUN 86 Previous editions are obsolete —— = SEGAL ® Eh AS SIG 
S/N 0102-LF-014-6603 Unclassified 


ae 


Approved for public release; distribution is unlimited. 


Why Stochastic Modeling is Essential in Analyzing 
Interdicted Transportation Network Performance 


Choi, Byung Kwon 
Major, Republic of Korea Army 
B.S., Republic of Korea Military Academy, 1980 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MANAGEMENT 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1989 


| P m 


ABSTRACT 


In this thesis, we use simulation to study the probabilistic dynamics of a road 
transportation network when the network is subjected to interdiction by aircraft. 
We consider several road interdiction schemes. The simulation results are compared 
to those obtained by using deterministic expected value optimization techniques. 
This latter approach has been used by other researchers in studies of interdictions of 
networks. We demonstrate, through the use of two realistic military examples, that 
the deterministic results poorly predict the performance of the stochastic systems 
and that the bias incurred by using deterministic methods is significant. Therefore, 


the stochastic model should be used in the real situation. 


lil 


TABLE OF CONTENTS 


INTRODUCTION se risis ooh prec waa = oe ce So st Neeser 1 
LITERATURE REVIEW 22 oe >> > oro cca eee een 4 
PROBLEM FORMULA TION =~ csceisnciete eicicisis > = eee 8 

A. NETWORK DESCRIPTION <2 cet cee > «+ eke g 
B. DETERMINISTIC PROBLEM FORMULATION. -:----+---++++++-- 11 
C. SOLUTION TECHNIQUE: ---------2-es eee e eee e cere eee ee eeees 11 
D. EXAMPLE «+ + cee oo oatears chee oo = 6 ete 6 13 
STOCHASTIC NETWORK MODEL---+-*>---+--- >>> >> - 16 
DISCUSSION OF SIMULATION PROGRAM: --+----+++eeeseeeeeeeee 19 
A. PROGRAM OVERVIEW: <> 2eeess 2-00 00s ee +056 - os en 19 
B NETWORK DATA FILPeces <7 meee ees oes - eos ee 90 
C. RANDOM NUMBER GENERATOR.----.----------+-0-++ =e 90 
D. GENERATION OF A RANDOM PROBLEM::--:---++--++++-+++++: J) 
FE. STATISTICAL PROCEDURES 222 oe ee eee 99 
TWO EXAMPLE- << «« 00 2 oe cere een tinct eee 93 
A. FEATURES COMMON TO BOTH EXAMPLES: -------+-+++-+:- 23 
B. EXAMPLE 1 (KOREA)----+++:eeceece eter teeter eee e teen ee eeces 95 
]. Recovery Times ccc ccc ccccescserccscsesccccrcccsceccccscccs oi 

2. Transit Time and Cargo Requirement----------+++++++++++-- 28 

3. Transit Time Comparison Between Two Models-----------+-+:: 29 


lV 


VII 


A. 
B. 


EXAMPLE 2 (CALIFORNIA): -++ee+seec eee ee rece eee e reece ee eee 30) 


1. Recovery Time---------+e eee e cece cece cece cece cece eenee 30 
9 SBpaweh TMP Ee 6 ocuuuacscl el OO eae 39 
3. Transit Time Comparison Between Two Models--------------: 32 
UEP NE By E2/2\ [SLC TY Oey Oy Re ERNE ok Ce ai 39 
Si NEN aIN Da GONGMUSIONG: ©. «ieee eaiiwwesces cares oe. 3A 
SUM AIPA 2 o 2 edo 6 06 DRI ee nn 34 
OUST CTLICE SEQUINS 6 Site 6 0 eC) nee en eae oe 35 


APPENDIX A: RANDOM CAPACITY IN MAXIMUM FLOW PROGRAM.: 37 


APPENDIX B: STOCHASTIC SIMULATION COMPUTER PROGRAM.-.--- 40 


APPENDIX C: TABLES AND NETWORK DATA FILE (EXAMPLE 1)--- 48 


APPENDIX D: TABLES AND NETWORK DATA FILE (EXAMPLE 2)--- 54 
SPIO ET ISRS ROE NGENG.. -. csc cs ccc ec cnc cc cc eee cee ete ce ete enee 58 


PEI CHIRISURION [IGT 020-6. - 2 ssc cs ec ccc ce eee ce cee ee eee ee 60 


LIST OF TABLES 


1. RANDOM NUMBER INTERVALS: «ec eee creer cece reece reese eeeeees p.)) 
2. LEVEL OF DAMAGE AND REQUIRED TIME+:++++e+eeese seer eee: 24 
3. OUTPUT OF RANDOM MAXIMUM FLOWS.-++++eesseeee rere eeres 38 
4. RECOVERY TIME DATA (KOREA) eee eee oe oo ier 8 48 
5. TRANSIT TIME DATA (KOREA)-++reeete rete etree teeter e cent cee: AQ 
6. NETWORK DATA FILE-<«-cceee eter ee cee tees eer errr ter teeeeeeerees 50 
7. TRANSIT TIME COMPARISON BETWEEN TWO MODELS 
(SCENARIO L)ss scree eee e cece ence net e nee teen nee e eens 52 
8. TRANSIT TIME COMPARISON BETWEEN TWO MODELS 
(SCENARIO 2)-++cceee cece cece een cece n eee een eet e eee Bo. 
9. RECOVERY TIME DATA (CALIFORNIA)----+ ++ sees sete este ee eee 54 
10. TRANSIT TIME DATA (CALIFORNIA): <*>: ++ s0++++ see) > oem 54 
11. NETWORK DA T AMD atts 4 veci reuse asus > . oer as 
12. TRANSIT TIME COMPARISON BETWEEN TWO MODELS 
(SCENARIO 3)--:escec secre cece eee ene n ee nent nnn ees 56 
13. TRANSIT TIME COMPARISON BETWEEN TWO MODELS 
(SCENARIO 4). +e ec cece cence cece eee e eee ene e een ee es 57 
14. DECREASE RATE AND CONFIDENCE INTERVAL:-----++-+++-++ +>: 33 


Vv] 


-— i 


A 


LIST OF FIGURES 


An Example Network--++++--++-see esse secre eect eter eee cence ecees 13 
Time-Dependent Network-+-+++s+ree eer e eee e ee eee eee eee n eee nees 14 
Maximum Flow Example Network-----+--+-+++seeceee eect eee eee eeee 37 
Example Network (Korea) SON Gano. 10 oO One 0 0 0 CREED C6 Cece eee 26 
Example Network (California): +-++++++ssseese cece ence eee e cence nee. 3] 


Vil 


ACKNOWLEDGEMENTS 
I wish to express my gratitude and appreciation to Professor Mike Bailey and 
Professor Allan W. McMasters for the instruction, guidance and advice throughout 
this research. 
Finally, many thanks to my wife, Dong—Hae and my daughter, Jung—Min, 
for their love, patience, and support during the past one and half years in Monterey, 
California. 


This paper is dedicated to my daughter Jung—Min. 


Vill 


I. INTRODUCTION 


In conventional or limited warfare the success or failure of a military 
campaign is directly dependent on abilities to transport material some distance 
overland within the required delivery time. This is particularly true in the case of 
armed aggression by one nation upon a contiguous or near — contiguous nation, such 
as might happen in the Korean Peninsula. 

The goal of military road network managers is to minimize the time required 
to transport war materials to their front line troops. The opposing or interdicting 
forces try to interdict some of the roads to maximize the transit time between the 
supply source and the front lines by using a variety of means, particularly airstrikes. 
If the transit time is greater than the required delivery time, the front line troops 
must either curtail their activity or retreat from the frontal area. Accordingly, it is 
to the advantage of the interdicting forces to delay the movement of materials as 
long as possible. 

The interdiction forces' target planner knows that the time when a target is 
hit is as important as which target is hit. In addition, partial destruction of a target 
using fewer resources may be sufficient for the immediate objective instead of total 
destruction with an attendant larger number of resources. 

Because the application of tactical air interdiction is probably the most 
effective method of denying the enemy vital war supplies, we will consider airstrikes 
as our only interdiction means. 

The extent to which a target is destroyed is unknown when airstrikes are 


planned. Aircraft may be shot down by anti—aircraft fire of the road managing 


forces, they may be prevented from completely destroying the target by defensive 
gunnery Or missiles, or they may fail to completely destroy the target due to human 
error. Therefore, the level of the success of an attack is probabilistic. As a result of 
the probabilistic success of attack, repair times, bypass construction times, and 
transit times after interdiction will also be probabilistic. We will assume that , once 
the probability of success of attack is known, these three times will be known for 
certain. 

In this thesis, we study the probabilistic dynamics of a road network 
transportation problem through the use of simulation methodology. We consider 
several road interdiction schemes and compare these schemes both by using 
deterministic expected value optimization techniques and by simulation. We will 
demonstrate, through the use of two realistic military examples, that the 
deterministic results poorly predict network transit time of the stochastic network 
and that the bias incurred by using deterministic methods is significant. 

Chapter II briefly discusses the assumptions, objective functions, and 
algorithms of earlier interdiction techniques and models. 

Chapter III includes the network description, problem formulation of a 
deterministic model, the solution procedure for that model and an example. 

Chapter IV describes the stochastic network model for which the transit time 
between source and sink is dependent upon the probabilistic success of an 
interdiction attack. 

Chapter V explains details of the simulation program such as the random 


number generator and the generation of random transit times. 


Chapter VI presents two realistic networks as examples, explains the 
network data file for each in detail, compares the results from the deterministic and 
simulation method, and conducts a sensitivity analysis between two models. 


Chapter VII presents a summary and conclusions. 


Il. LITERATURE REVIEW 


Current literature in the field of networks contains several different measures 
of effectiveness for interdicting a military transportation system. ‘These measures 
include the maximum delay of material transit, the greatest reduction to the 
maximum flow, and the least accumulated flow over the specified operational 
period. All of these measures involve single source to single sink, single commodity 
flows. 

Nugent [Ref. 1] presented a method of solving the problem of allocation of 
effort in the interdiction of a transportation network under the assumption that the 
damage function was deterministic and exponential. For example, even if a bridge 
is destroyed, it may still be possible to ford or ferry supplies across a river. The 
exponential damage function exhibits a decreasing marginal return from interdiction 
efforts. In addition, he assumed a planar network, divisible effort, single time 
period, and no storage or depot capability. His objective function was to maximize 
the reduction in the enemy's resupply capability subject to a constrained number of 
simultaneous airstrikes. Nugent used the maximal—flow minimal cut theorem and 
the topological dual to determine which route to select. He did not consider the 
probability of success of an attack. He used the expected values of flow over a given 
road segment. As we shall indicate, this use of expected flows detracts from the 
accuracy Of Nugent's results and methods. Finally, he assumed the road managing 
force would repair the interdicted roads so that they are operational some time after 


an interdiction. 


Sullivan {Ref. 2] developed a method to maximize the time required for an 
enemy to deliver war material to the front line troops. Any attack was considered 
completely successful, i.e., the target was fully destroyed. Repair and bypass 
construction times were assumed to be known linearly decreasing functions of time. 
He assumed that the resulting transit time over repaired arcs or constructed 
bypasses was always more than the uninterdicted transit times. 

Assuming certain success of all interdicting actions, he was able to solve a 
deterministic time—dependent shortest route problem using the method of Cooke 
and Halsey [Ref. 3]. The solution method used by Sullivan was to select a set of 
arcs to attack based on marginal increases in the source-to—sink delivery delay 
time. 

Mustin [Ref. 4] assumed that upper and lower limits on road capacities, as 
well as the amount of reduction per sortie, were known deterministically, and the 
reduction in capacity per sortie was linear between the upper and lower limits of 
capacity with the latter being positive. The goal of interdicting forces was to 
prevent the build—up of supplies over a 24—hour period. A computational procedure 
was developed for determining the optimum strike plan for minimizing network flow 
capacity. The work of Mustin, Sullivan, and Nugent were all masters’ theses at the 
Naval Postgraduate School under the direction of Alan McMasters. 

Wollmer (Ref. 5] presented two algorithms for targeting strikes in a 
lines—of—communication (LOC) network. He assumed that the user of the LOCs is 
attempting to achieve a circulation flow at minimum cost, a very general goal that 
includes, as special cases, maximizing flow between two points, meeting required 
flows between two points at minimum cost, and combinations of these two. One 


possibility is to assume all strikes allocated against a particular arc are directed 


against the same target and each has an identical and independent probability of 
successfully destroying it. The expected result of the attack is then used in an 
optimization scheme to maximize the minimum cost of communication. 

His algorithms attempt to make such arc costs as large as possible over time 
while decreasing arc capacities. The first algorithm considers arc costs as linear 
functions of flow; the second considers arc costs as piecewise linear functions of flow 
with one break point. Allocation is done on the basis of immediate user cost, repair 
times, and repair cost. Specifically, if a single strike is to be targeted, it is directed 
against an arc of maximum strike value, where strike value is defined as the repair 
cost plus the resulting cost increase per unit time of a minimum-—cost circulation 
flow multiplied by the repair time. For multiple strikes, no method for allocating 
strikes optimally, other than complete enumeration, is known. Wollmer solved the 
problem heuristically by repeated application of the one—strike algorithm. 

In Sengoku [Ref. 6], two new indices measuring the degree of influence of the 
arc on the vulnerability of the network were studied. He also assumed the 
interdiction of each arc was deterministically successful. 

The Network Interdiction Model (NIM) [Ref. 7] is a strategic military 
decision support system which treats three criteria; the maximum delay the 
interdicted targets would cause, the greatest reduction to the maximum flow, and 
the least accumulated flow over the specified operational period. 

Most of the approaches discussed above were concerned with "deterministic 
expected values," that is, the expected level of success of an attack was designated 
deterministically. Some of them assumed that the targets are struck one at a time 
sequentially. This thesis will examine a stochastic system which incorporates the 


probability of success of an attack and simultaneous attacks. In the next chapter, 


the algorithm for determining the time-dependent lengths of arcs will be developed. 
This algorithm will be used later in both the deterministic and stochastic network 


models. 


I. PROBLEM FORMULATION 


A. NETWORK DESCRIPTION 

The transportation system can be represented by a network of arcs and 
nodes. Nodes represent interactions of road segments. Further, they may be used 
to represent any point at which it is convenient to distinguish between the road 
characteristics on either side of the node. Arcs represent road segments. Each arc 
joins two nodes; that is, each is characterized by a beginning node and an ending 
node. The network is assumed to have undirected capacitated arcs. 

Two special nodes are the source and sink. The network is assumed to have 
a single source through which flow enters and a single sink through which flow 
leaves. If there are several sources or sinks, this may be taken care of by adding 
artificial nodes and arcs. The capacities of these artificial arcs are infinite and they 
are not vulnerable to attack. 

We are interested in the transportation times between the source and the 
sink after several arcs in the network have been attacked and destroyed. 

We will follow the network descriptions given in the Sullivan thesis[Ref. 2]. 
The notation (i, j) represents the arc between nodes i and j. Nodes are numbered 
from 1 to n with 1 for the source and n for the sink. The intermediate nodes have 
any of the values between 1 and n. There exists a transit time, di for each arc 
based upon road conditions, terrain, and time of day. Engineering units are 
stationed at several prespecified nodes in the network. These units must first travel 


to an attacked arc, and then perform the repair or make a bypass. Thus, each arc 


will have attributes of time to move between the closest engineering unit location 
and the damaged arc, time to repair, or time to construct a bypass. Let 


m= time to move between the engineering unit location and the damaged 


ana (iran) 

ee time to repair the destroyed structure in arc (i, j), 

bi; = time to construct a bypass in the arc (i, j), 

2 eae [or, b: time to recover the destroyed arc (i, j), 


c(t) = time remaining to recover arc (i, j) at time t, 
di,(z) = transit time over arc (i, j) once it has been repaired, 
d. b>) = transit time over a newly constructed bypass of arc (i, j) 

The values of di (7) and d;,(b) are assumed to be greater than or equal to the 
uninterdicted transit times, qi: This is a reasonable assumption since bypass 
construction is usually inferior to the original road segment, thus slowing traffic. 
Furthermore, even with a segment fully repaired, a vehicle may have a tendency to 
travel at a slower pace in an area that has recently undergone a bombing attack. 

The transit time can vary with types of military column formations. <A 
fundamental for convoy command and control is that the column be organized to 
meet the mission requirements and provide the degree of control necessary. ‘The 
convoy commander decides how his column will be organized for control, choosing 
from three basic methods: Close Column, Open Column, Infiltration. 

The difference in formation depends largely on vehicle spacing. The number 
of vehicles (density) per kilometer of road and the rate of march are accepted 
numbers (values) for average conditions when a movement is not influenced by 


attacking forces' actions. However, weather, tactical situation, attacking forces' 


capability, condition and type of road, vehicular maintenance, types of vehicles, and 


command policies may cause changes in average densities and speeds.{Ref. 8: p. 54] 
This thesis will assume the open column which is applicable for a normal supply 
support situation. 

Although any segment of a road is subject to attack, a highway segment 
that, after attack by strike aircraft, will be extremely difficult to repair or bypass is 
called a choke point. In other words, choke points are segments in a network which 
once attacked, force an enemy to either reroute traffic or expend large amounts of 
resources to keep the attacked segment open to traffic. Typical choke points would 
be bridges, mountain roads, and tunnels. In case no choke points exist in an arc, 
normal roads are selected for the targets. Each arc is assumed to have a unique 
attractive target for interdiction. Bypass and recovery times are based on the 
extent of the destruction of this target. 

Immediately following interdiction of a road segment, the arc transit time 
value becomes a function of recovery time. Therefore, after interdiction, the transit 
time over arc (i, j) will be defined by: 

dii(t) Sine ic; (t) 42 di(r), Gj (t) 4° di,(b)] 
where t represents the earliest clocktime after attack that a vehicle can depart from 
the source and arrive at the point of interdiction. The function which is minimized 
guides the road network manager in deciding how to get the arc back into operation. 
If c.(t) + di(r) < c..(t) + d. 


] J 1) j 


road segment. Conversely, if c(t) + di,(t) > c(t) + di,(b), then the manager will 


choose to construct a bypass around the interdicted point. 


(b), then the road manager will choose to repair the 


An attack is not always successful. An attack may have no effect, it may 
leave the road damaged but passable, it may cause the road to be impassable with 


minor repairs, or it may completely destroy the target, necessitating lengthy repairs 


10 


to restore the road's usefulness. We will call these levels of attack success no 
destruction, ruts, partial, and complete, respectively. Accordingly, the recovery time 


and the transit time after recovery will be different. 


B. DETERMINISTIC PROBLEM FORMULATION 

Road network managers want to minimize the time required to transport war 
material to his front line troops while attacking forces try to maximize the time 
between source and sink by interdicting this supply operation with airstrikes. 

If we assume that the success of each airstrike is known, we could solve a 
deterministic optimization problem which would give us the length of the required 
transit time. One method for solving this problem is linear programming, a well 
known method for performing optimization on networks. However, there exists a 
special purpose algorithm for this particular optimization problem due to Cooke and 


Halsey {Ref. 3]. This method was used by Sullivan [Ref. 2]. 


C. SOLUTION TECHNIQUE 

To determine the shortest transit time for a given set of points of 
interdiction, the Cooke and Halsey algorithm compares the repair and construction 
bypass functions for each damaged segment along the route. The "length" of 
interdicted arc is then replaced with d; jt). The arcs are assumed to be attacked 
simultaneously within available air sorties. 

The algorithm is given in the following steps. 


1. Choose the arc sets to attack. 


1] 


2. The transit time of the arc to be interdicted is replaced by the time 


dependent length function 


aI) 5 


di(t) = c(t) + dj, 


j 
if repair is optimal; otherwise, 
di,(t) = c(t) + di,(b). 

The problem is now time dependent in nature. 

3. Define the new tentative node label f. to be an upper bound on the 
earliest time of arrival at node i, and permanent label, iP, to be the earliest possible 
(optimal) time of arrival. 

4. Permanently label node 1 (source) with a value of f, = 0 and label all 


other nodes with values of infinity; i. e., 


£ (0) oP, = =f) = 


l eee n OO. 


5. Tentatively label all nodes j with the minimum of the current node label 
f and the sum of f, and dij «(11.0% 
it ) = min [f, +d 1h PY. (0) 7. 


J 
6. Find the aa tentative node label; i.e., f,, and declare it permanent. 
f,P = = 0gu i (1 ) 
eal. 


7. Node k, the new permanent node, is then used to attempt to reduce the 
labels at all tentatively labeled nodes by comparing tie Sh qi (f, ) to the current 
label. The minimum new temporary label is declared permanent and used as a basis 


for the next iteration; i.e., 


fm) = min [f,P + 4,.(f,), pmo} |: 
{P= = hin fm m)_ 
j¢#k ? 


§ Terminate when node n is permanently labeled. 


If the network is composed of N arcs and s arcs can be interdicted, there 
exist . possible interdiction sets. The interdiction set which has the maximum 
value of ‘ae contains the optimal points of attack for the transportation network so 
as to maximize the time required for the road manager to transport war material to 
his front line troops. This is the optimal solution to the dual of the road manager's 


problem. 


D. EXAMPLE 
Consider the simplified uninterdicted transportation network described by 
Figure 1. All transit times ( di, ) are in minutes, and nodes 1 and 6 will be the 


source and sink nodes, respectively. 





60 


Figure 1. An Example Network. 
We believe that the most realistic planning model involves the assignment of 


a number of aircraft sorties to different network arcs so as to perform a single, 
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simultaneous airstrike. Thus, for our example it is assumed that exactly two arcs 
can be attacked simultaneously at time 0 with allocated air sorties. 
1. The arc set, (2, 4) and (3, 5), will be interdicted. 
2. f,() = (1) — 30, f,(?) =o) 
3. The arc to be interdicted is replaced by the time function. Figure 2 shows 
the results. 
Suppose that arc (2, 4) has the following data: 
Mo, = 20, 
fo, = 180, do, (r) = 60, 
bo4 = 200, do,(b) = 70. 
Next, Mo4 + To4 + do4(t) = 260 < mo4 + bog + do4(b) = 290. 
50, Cog = Mog + To4 = 200. ( the road manager will choose repair. ) 






d,,(") = 60 


I5 =280 


Figure 2. Time—Dependent Network. 
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According to the time-dependent function, 


oy yeoman ch iow in liind., (a= 230: 


Next, suppose the data for arc (3, 5) is 


Mg, = 30, 
ie = 200, dg-(r) =a 
be. = 180, dg-(b) 5), 


Then, mg, + Ign + dg-(r) = 300 > Mgn + ben + dgn(b) = 290. 
50, Con = Man + bg. =i: 
Car(f3) =O) do-(f3) = 2. 

4. By using the above results, £,{3) = 260, £,(4) = 280, £5) = 290. 

The road manager would send his supplies over the shortest route after 
interdiction. For this example that route is 1 — 2 — 4 —6. Before sending these 
supplies he must repair arc (2, 4). He can send supplies out from node 1 at time 
zero but the trucks will have to wait at node 2 for 170 time units until the arc is 
repaired. The other alternative is to wait at node 1 for 170 time units before 
departing to node 2. If this tactic is used then the arc (2, 4) will just be repaired as 
the trucks arrive at node 2. 

In the next chapter, we will consider a model in which the success of a given 


attack is random, and observe the effects of these attacks on source—to—sink transit 


time. 
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IV. STOCHASTIC NETWORK MODEL 


The assumption that any attack is considered completely successful is far 
from the real situation. The target damage may be complete, partial, or only ruts 
in the road may result from a successful weapon delivery. Also aircraft may be shot 
down by ground anti—aircraft gunfire or by defensive air patrol interceptors of the 
road managing forces. Probabilities of kill on the target vary with the levels of 
pilots' skill, the opponent's defensive gunfire, the exposure time of an aircraft, the 
type of munition, hardness of target, and so on. Therefore, the level of the success 
of an attack must be modeled as probabilistic or random instead of assuming that 
we always have completely successful destruction of the target. 

As a result of random successes, repair times, bypass construction time, and 
transit times after interdiction will also be random variables. Repair times and 
bypass times will be designated according to the probability of success of attack, 
i.e., the level of damage on the target. In case of no damage or ruts only, the 
construction time will be 0. The values after complete destruction are assumed to 
be greater than or equal to the values after partial destruction. Also, the transit 
time after interdiction will be influenced by the success of the attack. Even though 
the interdicted arc is recovered fully, a convoy column commander may have a 
tendency to travel at a slower pace in an area that has recently undergone a 
devastating bombing attack because of the psychological effects of being 
demonstrably vulnerable. It is assumed that the transit times associated with the 


rut case are larger than the other transit times. This is a reasonable assumption 
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since the road manager will probably not repair a roadway having only minor 
destruction. 

To reiterate, the transit and recovery times of the attacked arcs are random 
variables which are a function of the effectiveness of the attack. Once the attack's 
success is known, the transit time and the reconstruction time will be known with 
certainty. 

When faced with random network arc performance, most analysts decide to 
take the expected value of the arc performance and use this expected value as if it 
were a deterministic quantity. The expected value is thus input into a deterministic 
network optimization scheme such as the Cooke and Halsey algorithm above, and 
the result is reported as the expected value solution. Several of the papers cited in 
the literature review, particularly Wollmer [Ref. 5], have such published results. 

Using expected values in cases such as maximum flow or shortest path 
problems can be shown to give extremely inaccurate results. To show this we will 
solve a network maximal flow problem using the expected values and the simulated 
results. The details are presented in Appendix A. For the simulations we used 
coins and dice to generate random capacities and calculated the resulting random 
maximum flows. We also calculated the expected value of the capacity of each arc's 
capacity and computed the maximum flow on the network using those expected 
capacities. 

As we can see from Table 3 of Appendix A, the result of using the expected 
Capacities is to produce incorrect solutions. The amount of flow calculated using the 
expected values is 10 and is greater than the result of the average of 6.27 from the 


stochastic simulation after 30 experiments or trials by a factor of 37.3 per cent. 
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Using expected values in place of random variables would lead the transportation 
manager to overestimate his ability to transport goods. 

At this point, we can predict similar results for the transit time. We expect 
that the transportation forces will be able to transport goods to the front much 
faster than we would be lead to believe had we used the expected time of transit for 
each arc. This result implies that the algorithms used in all of the deterministic 
models described in our literature search may not be as efficient as reported. The 
goal of the remainder of this thesis is to demonstrate this for the Cooke and Halsey 
algorithm. Chapter VI will demonstrate that the expected transit times are bigger 


than values calculated using the simulation program we will describe in Chapter V. 
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V. DISCUSSION OF SIMULATION PROGRAM 


If the relationships which comprise the model are simple enough, it may be 
possible to use analytical methods (such as algebra, calculus or probability theory) 
to obtain exact information on questions of interest. However, most real—world 
systems are too complex to allow the models to be evaluated analytically, and these 
models must be studied by means of computer simulation. In a simulation we use a 
computer to evaluate a model numerically over a time period of interest, and data 
are gathered to estimate the desired true characteristics of the model. [Ref. 9: p.1] 

A simulation model is said to be deterministic if it contains no random 
variables. For a deterministic model, there is a unique set of model output data for 
a given set of inputs. On the other hand, a simulation model is stochastic if it 
contains one or more random variables. The output data for a stochastic model are 
themselves random and thus only estimates of the true characteristics of the model. 
(Ref. 9: p.3]. Our model will be stochastic because the attack successes are random. 
The output data we will collect are the transit times from the source to the sink. 
thus our output is indeed random. We will compute statistics based on this data 
such as the sample mean and sample variance. Our program allows us to specify the 


number of statistics we desire. 


A. PROGRAM OVERVIEW 
We wrote our simulation program in Professional FORTRAN, by Ryan 
McFarland. The program was implemented on an IBM PC compatible computer. 


Our program begins by inputting the network data file which contains all of the 
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relevant information for the problem. We then enter a loop in which the 
pseudorandom attack successes are calculated. The resulting deterministic 
time—dependent shortest path problem is generated and solved using the Cooke and 
Halsey algorithm, and the resulting shortest path length is recorded. 

This loop is repeated the specified number of times. From the recorded 
statistics, the sample mean and sample variance are calculated and an output file is 
built which contains the problem, iteration results, and the statistical results. 
Appendix B contains a listing of the simulation program. We now provide some 


more detail concerning each stage of the program. 


B. NETWORK DATA FILE 

As we can see in the Appendix C, the network data file (Table 6) includes 
the following information concerning the network itself: source, sink, number of 
nodes, number of arcs, the tail (T) and head (H) node of each arc. For each set of 
attacked arcs, we specify the recovery times for each type of destruction for each 
arc (c(1), c(2), c(3), c(4)), the transit time before destruction (d(0)), the transit 
times after recovery according to the success of attack (d(1), d(2), d(3), d(4)), the 
arc set attacked, and the probability of each type of destruction is given for each arc 
(p(1), p(2), p(3), p(4)). Thus, for each attacked set, we have the option of changing 


any relevant recovery or transit time data. 


C. RANDOM NUMBER GENERATOR 
Because the success of attack is probabilistic, we utilize pseudorandom 


numbers, which is conceptually equivalent to flipping coins or rolling dice. This 
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stochastic model will use the linear congruential method (LCM) which is one of the 
most effective and popular types of random number generators. The FORTRAN 
function RAND() of the program is for generating random numbers. 

LAST is the seed which is a locally saved variable. Its value is updated 
throughout execution of the subroutine that uses RAND. LAST has the value that 
depends on its previous value and the formula involving LAST, C, and xM. The 
value of LAST will be between 0 and xM-—1 (inclusive) and will be suitably random 
for proper choices of the constant values xM, C and the starting value of LAST. 
The values of 29, 217, and 1024 for LAST, C, and xM, respectively, as in RAND, 
define a usable LCM random number generator. |Ref. 11: pp. 175-178] 

With these values of LAST, C, and xM, the function RAND will generate 
pseudorandom numbers which are approximately uniform on the interval {0, 1). 
Because of our relatively small sample sizes, cycle length was not an issue. As an 
experiment, we also used the recommended generator in Press, Flannery, Teukolsky, 
and Ve Herling [Ref. 11]. Our results did not differ significantly using these values 
for LAST, C, and xM. 


D. GENERATION OF A RANDOM PROBLEM 

We will use different values for different arcs because each arc has different 
distance, defensive posture, hardness of target, and so on. From these known 
probabilities we can then construct the associated random number intervals. This 
process is Summarized in Table 1. Suppose that the random number 0.47 is 


generated, then Table 1 would give the result that the arc was partially damaged. 
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TABLE 1. RANDOM NUMBER INTERVALS 


Level of Probability Cumulative Interval of 
Damage (each level) Probability random number 
0 0.00 — 0.20 


0.21 — 0.40 
0.41 — 0.90 
0.91 — 1.00 


If the arc sets to attack are decided upon, the probability of success of attack 





will be designated from using the random number generator and Table 1. 
Accordingly, the data file then provides the recovery time and the transit time after 
recovery. Finally, by using the Cooke and Halsey algorithm, we can get the transit 


time from source to sink. 


E. STATISTICAL PROCEDURES 

At the end of each series of simulation runs, the subroutine STATS 
calculates the sample mean and sample variance as well as the 90% and 95% normal 
confidence intervals for the source—sink transit time. Note that this implies that 
the simulation model should not be used with less than 30 iterations if one truly 
desires normal confidence interval results. We determined that, with 200 samples, 
most of our confidence intervals for the mean transportation time were reasonably 
tight. Although increasing sample size would have tightened these intervals further, 


we chose to keep the sample size the same for all cases in the interest of uniformity. 
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VI. TWO EXAMPLES 


In this chapter, we will compare the deterministic model and the stochastic 
model for two realistic situations: the road network between Kaeseong and Sariwon 
in North Korea (Example 1) and the road network between Gilroy and Carmel 
Valley Village in Central California (Example 2). We made reasonable assumptions 
about moving times, bypass and repair times, and transit times as input into the 
simulation program. This input we call the network data file. We will first explain 
the procedure of constructing the network data file. We will then compare the 
outputs of two models, analyze the results, and interpret what the differences 


between two models represent. 


i FEATURES COMMON TO BOTH EXAMPLES 

It is assumed that the probability of each level of destruction is known and is 
the same for each arc. The values we selected were shown in Table 1 of the last 
chapter. 

It is assumed that the size of the operational unit in both examples is one 
North Korean infantry division. The theater—level average consumption rates in a 
temperate zone are given as 123.7 pounds per person per day [Ref. 12: p.2—5]. If we 
assume that the number of personnel of the North Korean infantry division is 
10,000, the amount of cargo needed to support such a division would be 1,237,000 
pounds per day. 

Let's assume that the North Korean Army uses trucks which have the same 


payload and cubic capacity as the South Korean cargo truck, M35A2, 2 1/2-ton. 
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Its payload is 10,000 pounds on a paved highway and 5,000 pounds offroad [Ref. 8: 
p.1—7|. We will consider only the highway case in this example. The number of 
trucks required to move 1,237,000 pounds is 128. 

We also assume that the transit time once an arc has been repaired and the 
transit time over a newly constructed bypass are equal, thus d, j0>) oats jo for each 
arc (i, j). Although the level of damage, hence the recovery time will be different 
for each arc, we assumed damage levels which were the same for each arc in our 
examples. In any operational use of this methodology, we must account for varying 
levels of damage. Table 2 shows the common degree of damage and recovery time 


in case of full destruction for each type of target. 


TABLE 2. DEGREE OF DAMAGE AND REQUIRED TIME 


Degree of Damage Required Time 
Bridge One pier and Panel Bridge : 5 hr 
two spans Bypass Construction : 3 hr 
Tunnel Demolition 50 m Repair : 20 hr 
from the entrance Bypass : impossible 


Mountainous Demolition 100 m Repair : 6 hr 
Road Bypass : 7 hr 
Road Demolition 150 m Repair : 5 hr 

Bypass : 4 hr 


Network data files, Table 6 in Appendix C and Table 11 in Appendix D, 





include the multipliers for inflating the transit and recovery times. For example, 
the recovery time will be assumed to be zero (c({1) = c(2) = 0) in case of no damage 
or ruts only because the arc is used without repair or bypass. We assume that the 
recovery time for partial destruction (c(3)) is 60 % of the recovery time for complete 


destruction (c(4)). We can get c(4) from Table 4 in Appendix C and Table 9 in 
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Appendix D. As for transit time, we get d(0) from Table 5 in Appendix C and 
Table 10 in Appendix D which is the transit time before interdiction (how these 
times were determined will be described below). Though there is no damage, the 
transit time will be increased due to psychological effects. In case of ruts only, we 
can expect considerable delay because road manager forces do not repair or 
construct a bypass. Furthermore, even with a segment fully repaired, a vehicle may 
have a tendency to travel at a slower pace in an area that has recently undergone a 
bombing attack. Consequently, we consider that d(1), d(2), d(3), and d(4) are 
increased by certain multipliers. (Refer to the end of Table 6 in Appendix C and 


Table 11 in Appendix D for details) 


By EXAMPLE 1 (KOREA) 

Current operational planning for the defense of South Korea is focused on the 
possibility of a North Korean surprise attack. It is estimated that North Korea may 
maintain assault pressure with the prepared amounts of supplies in the frontal areas 
for three days. On the fourth day, additional supplies would be needed for the 
continued assault. Let us assume that it is the fourth day of the assault. We will 
consider the network between Kaeseong and Sariwon to be as shown Figure 4. As 
we can see in Figure 4, the actual road network is composed of 26 nodes and 38 arcs. 
The source node is Sariwon and the sink node is Kaeseong. Kaeseong is near to the 
Demilitarized Zone (DMZ) and it is reasonably assumed that the North Korean 
logistics command is located in Sariwon. This area is mainly composed of 
mountains, rivers, and rice paddy fields. The road network is well developed in 


contrast to other areas of North Korea. 
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Figure 4. Example Network (Korea) 


Data concerning North Korean equipment for transportation and road repair 
is unavailable in-an unclassified form. In lieu of this data, we used the U. S. Army 
field manual [Ref. 8, 12. 13] to supply us with plausible data. We will now describe 
the development of the recovery and transit time data used in the network data file 
for the Korean scenario. 

1. Recovery Time 

The recovery time of an arc is the time required to repair a damaged 
road segment or to construct a bypass. We will assume that all network attacks are 
in the form of South Korean air strikes. We assume that eight sorties per day are 
available to attack the network and two sorties are allocated for each target. These 
are requirements which are imposed on strike aircraft operations in order to ensure 
safety of flight. Therefore, the number of arcs to be attacked is four. 

Let's assume that there are five engineering units in this area: Sariwon 
(node 1), Packchon—ni (node 9), Chongdan (node 16), Kumcheon (node 18), and 
Kaeseong (node 26). 

By using the assumptions stated above, the data in Table 4 of 
Appendix C was calculated. Column (1) of this table is the number of the arc (the 
nodes each arc connects to are shown in Columns 2 and 3 of Table 6). Column (2) 
represents the types of targets: bridges, mountain roads, tunnels, and roads. 
Column (3) is the distance in kilometers between the engineering unit and the 
target. Column (4) is the movement time calculated from column (3) and is stated 
in minutes. It is assumed that the speed of recovery personnel and equipment is 30 
KM per hour. Columns (5) and (6) are repair time and bypass times expressed in 


minutes. Column (7) gives the minimum of the previous two columns. Column (8) 


2/ 


is the recovery time which is the sum of columns (4) and (7). Column (9) gives the 
location of the nearest engineering unit to each target. 
Ze Transit Time and Cargo Requirements 

As we discussed in Chapter III, we will choose open column 
formations for traversing all arcs. The density per KM is 12 trucks and the speed of 
each truck is 24 KM/H [Ref. 8: p.5-5]. The additional assumptions are no 
exogenous traffic and concrete or bituminous pavement. 

Transit times are called road clearance times in military terminology. 
It is the total time a column or element thereof requires to travel over and clear 
either a section or all of a road. Road clearance time equals the time distance plus 
time length. Time distance is the time required to move from one point to another 
at a given rate of speed. Time length (pass time) is the time required for a column, 
or element thereof, to pass a given point [Ref. 8: p. F—2]. For application to the 
real network, the transit time over an arc is assumed to be only the time distance 
while the transit time of the last arc of each route, such as arcs 24 and 38, is the 
time length plus the time distance. This is due to the fact that the convoy column 
experiences the time length only at the end of their journey. Thus, only the last arc 
of each path needs to have time length added. The formula for time distance [Ref. 


8: p. F—9] is: 


Time Distance = Distance / Rate , 
The formula for time length [Ref. 13: p. 3-57] is: 
Time Length = (number of vehicle * 60) / (density * rate) + time gaps +EXTAL, 


where EX TAL is the extra time allowance. 


In the above formulas, we know all of the required data except for the time gaps. 
For time gaps, we need the vehicle gaps. In the open column, the density is 12 
trucks per KM and the length of a truck is 17 feet (1 feet is 0.3048 meters). Thus, 
the vehicle gap is: 
Vehicle Gap = (1000-12 * 17 * 0.3048) / 11 = 85 (M). 
In the open column, the traveling rate is 24 KM / H. So, the time gaps value is: 
Time Gaps = (127 * 85 * 60) / 24000 = 27 (MIN). 
The time length value is: 
Time Length = (128 * 60) / (12 * 24) + 27 + 4 = 57 (MIN). 
Table 5 includes the above information. Column (2) contains the length (Km) of 
each arc; Column (3) lists time distance; Column (4) lists time length; and Column 
(5) gives the transit time before interdiction. 
a Transit ‘Time Comparison Between ‘Two Models 
Given the arc sets to be attacked, we will compare the results of two 
models. The deterministic model used expected values of recovery times and transit 
times, 1.€., 
d(e) = 0.2d(1) + 0.2d(2) + 0.5d(3) + 0.1d(4), 
c(e) = 0. 


Outputs of the stochastic model was from 200 simulations. Table 7 of Appendix C 


2d(1) + 0.2d(2) + 0.5d(3) + 0.1d(4). 


shows the results of 10 arc sets selected as reasonable interdiction possibilities from 
the ri choices(scenario 1). As we suggested in Chapter IV, the average network 
transit times from the stochastic model were less than those for the deterministic 
model except for the first, second, and tenth arc sets. The main reason for the 
reversal for these three arc sets is that the recovery times have less influence on the 


time—dependent transit time when the interdicted arcs are near the sink. So, in an 
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attempt to increase the effects from such arcs, we multiplied recovery times by two 
and reran the simulation. The results are presented in table 8 (scenario 2) of 
Appendix C. As the reader can see in this table, the second and fourth arc sets still 
behave contrary to our expectations. The reason is that the arc sets include arcs 24 
and 38 which are the last arcs of each route. Recovery times have no influence on 
the time—dependent transit time of these last arcs because they are repaired before 


the convoy reaches them. 


C. EXAMPLE 2 (CALIFORNIA) 

The assumptions of Example 2 are similar to those of North Korea. 
Therefore, we will use most of the same assumptions, and same procedures. 

The area is between Gilroy and Carmel Valley Village. It is assumed that 
Carmel Valley Village is the frontal area and the logistics command is located at 
Gilroy. As we can see the Figure 5, the network is composed of 11 nodes and 15 
arcs. The source node is Gilroy and the sink node is Carmel Valley Village. 

ie Recovery Time 

The structures which can be choke points are only bridges and roads 
in this area. It is assumed that six sorties per day are available in this area and two 
sorties are allocated for one target. Therefore, the number of arcs to be attacked is 
three in any one arc set. 

Let's assume that there are three engineering units in this area: Gilroy 
(node 1), Castroville (node 6), and Carmel Valley Village (node 11). The data of 
Table 2 is assumed applicable in this area also. Table 9 of Appendix D shows 


recovery time data. The arc distances in Column (3) are expressed in miles. It is 


assumed that the speed of recovery personnel and equipment is 30 mile per hour 


(mph). 
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Figure 5. Example Network (California) 
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94, Transit Time 
For the open column, the density per mile is assumed to be 20 and the 
speed is assumed to be 15 mph [Ref. 8: p.5—5]. Table 10 of Appendix D shows the 
transit time data. Table 11 of Appendix D is the network data file which applies for 
this example. 
SE Transit Time Comparison Between Two Models 
Table 12 of Appendix D shows of 10 arc sets selected as reasonable 
interdiction possibilities from the (3° choices (scenario 3). The average network 
transit times from the stochastic model were less than those for the deterministic 
model except for the first and ninth arc sets. The main reason is the same to 
scenario 1. Table 13 of Appendix D is for the case recovery time data is double that 
of scenario 3 (scenario 4). In case of scenario 4, every selected arc set has the same 
tendency, the transit times from the deterministic model is more than those from 


the stochastic model. 


D. GENERAL COMMENTS 

In all of the scenarios we see the general tendency is that the expected 
network transit time of the deterministic model is more than the average times of 
the stochastic model except when the arcs to be attacked are near the sink. This 
latter result 1s because the accumulated transit time to a node for such arcs is 
usually more than the recovery time of the arc. Only when the recovery times 
become large do those arcs recovery times affect the transit time through the net 
work as we Saw in scenarios 2 and 4. We also notice that, in those cases where the 
deterministic expected time is larger than the average transit time from the 


simulations, the deterministic times are also larger than the upper bound of the 95 
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% confidence interval for the simulated times. Finally, we notice that the 
differences between transit times of two models is significant and the results are, at 
times, unpredictable. 

The average decreases in transit time from the deterministic model to the 
stochastic model and the width of the confidence interval are shown in Table 14. 

As we can see from these examples, the transit time calculated using the 
expected values is usually larger than the averages resulting from the stochastic 
simulation. Using expected values in place of random variables would lead the 
transportation manager to underestimate his ability to transport goods within a 
required delivery time. On the other hand, using the expected value model to 
determine the effectiveness of airstrikes to interdict flow of supplies leads to 
overestimating this effectiveness. In cases where this effective interdiction is 
critical, such as day four of a Korean conflict, this overestimation of effectiveness 


could cause serious losses of troops and territory. 


TABLE 14. DECREASE RATE AND CONFIDENCE INTERVAL 


Scenario | Avg. Decrease rate (%) Range Of Width Of C.I (95 %) 


3.1625, 37.4108 
3.1625, 17.8246 
1.7914, 18.6705 
1.7914, 36.6148 





VII. SUMMARY AND CONCLUSIONS 


A. SUMMARY 

The goal of military road network managers is to minimize the time required 
to transport war materials to their front line troops. The opposing or interdicting 
forces try to interdict some of the roads to maximize the transit time between the 
supply source and the front lines by using a variety of means, particularly airstrikes. 
If the transit time is greater than the required delivery time, the front line troops 
must either curtail their activity or retreat from the frontal area. Accordingly, it is 
to the advantage of the interdicting forces to delay the movement of materials as 
long as possible. 

The assumption that any attack is considered completely successful is far 
from the real situation. The target damage may be complete, partial, or only ruts. 
Also aircraft may be shot down by anti—aircraft gunfire including defensive air 
patrol interceptors of the road managing forces. Probabilities of kill on the target 
vary with the levels of pilots' skill, the opponent's anti—air capability, the exposure 
time of aircraft, the type of munitions, hardness of targets, and so on. Therefore, 
the level of the success of an attack must be modeled as probabilistic instead of 
assuming that we always have completely successful destruction of the target. 

If the relationships which comprise the model are simple enough, it may be 
possible to use analytical methods (such as algebra, calculus or probability theory) 
to obtain a direct optimal solution to questions of interest. However, most 
real—world systems are too complex to allow the models to be evaluated 


analytically, and the only resort left is to use computer simulation. In such a 
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simulation we simulate a single occurrence of the events of the situation being 
modeled. Repeated simulations provide data which can be used to estimate the 
desired true characteristics of the model. 

Our computer simulation program in Appendix B began by inputting the 
network data file which contained all of the relevant information for the problems of 
interest. We then entered a loop in which the pseudorandom attack successes were 
calculated. The resulting deterministic time—dependent shortest path problem was 
generated and solved using the Cooke and Halsey algorithm, and the resulting 
shortest path length was recorded. This loop was repeated 200 times for a given set 
of interdicted arcs. From the recorded statistics, the sample mean and sample 
variance of the resulting transit time was calculated and an output file was built 


which contained the results of the iterations, and the statistical summary. 


1. CONCLUSIONS 

The result of using the expected arc transit time values is to produce biased 
solutions. The transit time calculated using the expected values tends to be larger 
than the average values obtained from stochastic simulation. Using expected values 
in place of random variables would lead the transportation manager to 
underestimate his ability to transport goods within a required delivery time. On the 
other hand, using the expected value model to determine the effectiveness of 
airstrikes to interdict flow of supplies leads to overestimating this effectiveness when 


transit time over the shortest route is the main concern. In cases where this 


effective interdiction is critical, this overestimation of effectiveness could cause 
serious losses of troops and territory for either side. Therefore, the stochastic model 


should be used in the real situation. 


APPENDIX A. RANDOM CAPACITY IN MAXIMUM FLOW PROBLEM 


A. PROBLEM DESCRIPTION 
We performed a simple experiment involving a small maximum flow problem 
on the network shown in Figure 3. We used coins and dice to generate random 


capacities and calculated the resulting random maximum flows. 


9 C 4 


Figure 3. Maximum Flow Example Network 
For random capacities, we used a dime, a quarter, a penny, a nickel, and a 
die simultaneously to generate the values of a, b, c, d, and e,respectively. If a coin 
flip results in a head, the capacity is 10, 25, 2, and 5 for arcs a, b, c, and 
d,respectively. If a result is a tail, the capacity is 1. As for arc e, the value is the 
face of die multiplied by two. 
The expected value of the capacities are 5.5 for arc a, 13 for arc b, 1.5 for arc 


c, 3 for arc d, and 7 for arce. 
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B. PROCEDURE AND OUTPUT COMPARISON 

The maximum flow on the network using the expected capacities is 10. 
Table 3 shows the output of random maximum flows. z, is the single commodity 
maximum flow, and fi. is the i sample mean after i independent experiments or 


trials. Thus, 


] 
=() z,. 


TABLE 3. OUTPUT OF RANDOM MAXIMUM FLOWS 
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As we see from the table, fl, stabilizes somewhere near 6.25. The amount of 
flow calculated using the expected values is greater than the result of the thirty 
stochastic simulations by a factor of 37.3 per cent. Using expected values in place of 
random variables would lead the transportation manager to overestimate his ability 


to transport goods. 
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APPENDIX B. STOCHASTIC SIMULATION COMPUTER PROGRAM 


program thesis 


C 
c Parameter statement must be identical in each subroutine except 
C for RAND and STATS 
parameter(maxnod = 40, maxarc = 60, maxtyp = 5) 
C 
parameter(kseed = 3645) 
C 
real ctabl(maxnod, maxnod), dtabl(maxnod, maxnod), 
& c(maxarc, maxtyp), d(maxarc, O:maxtyp), cdf(maxarc, maxtyp), 
& ztabl(maxnod, maxnod), sum, sumsy, label(maxnod) 
C 
integer numnod, source, numtyp, numarc, atakd(maxarc), 
& head(maxarc), tail(maxarc), dest, kk, iter 
C 
character*20 fname 
characler*t#l 2 
C 
common last 
C 
C Initialize the seed, sums, and the data file name 
last = kseed 
fname = 'netwrk2.dat' 
sum = 0.0 
sumsq = 0.0 
do 1019 i = 1, numnod 


do 1019 j = 1, numnod 
ztabl(i1. 7) ~-=-0.0 
1019 continue 


C 
C Determine the number of simulation iterations to run 
C 
Print *, ‘How many iterations would you like, sir?’ 
print *#, ‘Iterations = ' 
read *, iter 
print * 
print * 
C 
Ci Input the network structure in the subprogram READNET 
C 
call READNET(fname, head, tail, c, d, numtyp, nunmnod, 
& numarc, atakd, cdf, source, dest) 
c 
Cc Debugging prints 
C 
G print * 
C print *, ‘Going from ', source, ' to ', dest 
C read (*#,1) 2 


do 1011 kl = 1, numarc 
if (atakd(ki).eq.1) then 
print *, ‘Arc attacked: tail = ', tail(k1), °' head = ', 
& head(k1), *' arc # ', kl 
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endif 
continue 


9191 is the big loop for the iterations 


do 9191 kk = 1, iter 


Generate random construction and travel time tables 
call RTABLE(ctabl, dtabl, tail, head, c, d, atakd, 
cdf, numtyp, numnod, numarc) 


Debugging print stalements 
print *, 'Dtable and Ctabl from the main after RIABLE’ 
print * 
do 1009 i = 1, numnod 
write(*, 10) (dtabl(i, j), J = 1, numnod) 
continue 
write(*, *) 
do 1008 i = 1, numnod 
write(*, 10) (ctabl(i, j), j = 1, numnod) 
continue 
read (*,1) gz 
format(al) 
Formaci10(fe.3, 2x), 7 1018.3, 2x) / 100f8.3,. 2x) / 
POt1TS8. 3, 2x)) 
print *, ‘The cdf for °', numtyp, ’types of destruction 
print * 
do 1010 i = 1, numarc 
print, (cdf(i, 3), jJ = 1, numtyp) 
continue 


Call the solver with the generated table 
call SPSOLVE(numnod, dtabl, ctabl, source, label) 
Debugging print statements 


print *, ‘destck ', dest, ‘here come labels’ 
do 23 kl = 1, numnod 
print *, label(k1l) 
continue 
print * 


Debugging print statements 
call outpt(label, numnod) 


Collecting statistics on the performance 


sum = sum + label(dest) 
sumsq = sumsq + (label(dest)**2) 
print *, ’Distance to source = ', label(dest) 


print * 


4] 


c print * 
Cc read (*, 1) gz 
9191 continue 


€ 
c Call the statistics subprogram 
C 
call STATS(sum, sumsq, iter) 
C 
stop 
end 
C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C 
subroutine STATS (sum, sumsq, iter) 


c 
real sum, sumsq, riter 
integer iter 
C 
riter = real(iter) 
print *, “your statistics sir... 
xbar = sum / riter 
print *, 'average = ', xbar 
if (iter.eq.1) then 
print *, ‘thats all we get for one iteration 
return 
endif 
stddev = sqrt((sumsq - ((sum ** 2) * (l/riter)))/(riter - 1)) 
print *, 'standard deviation’, stddev 
print *, ‘stdev of x-bar = ', stddev / sqrt(riter) 
print *, ’ for iter = ", iver 
radl = 1.96 * stddev / sqrt(riter) 
rad2 = 1.645 * stddev / sqrt(riter) 
print *, °'95% confidence interval (’,xbar - radl,’,’,xbar 
& + radi,')’ 
print *, '90% confidence interval (',xbar - rad2,',',xbar 
& f radZ, )° 
print * 
return 
end 
Cc 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCUCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C 
subroutine SPSOLVE(numnod, dtabl, ctabl, source, label) 


C 
parameter(maxnod = 40, maxarc = 60, maxtyp = 5) 

C 
integer numnod, source, curnod, node, prm(maxnod) 
real dtabl(maxnod, maxnod), ctabl(maxnod, maxnod), 

& curtim, label (maxnod) 

g 

c print *, 'Dtable and Ctable inside solver' 

C do 1009 i = 1, mnumnod 
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C write(*, 10) ldPebil (i, in) ee J = i numnod ) 
C009 continue 


c write(*, *) 
c do 1008 i = 1, numnod 
C write(#*, 10) (ctabl(i, J), J = 1, numnod ) 
C008 continue 
C read (¥,1) 2z 
1 format(al ) 
10 formate 1s.3, 2x), / 10(f8.3, 2x) / LO(fF8.3, 2x) / 
& 10(f8.3, 2x)) 
do 1000 i = 1, numnod 

prm(i) = 0 

label(i) = 99999.0 
1000 continue 
C 
c print *, 'SOURCE = ', source 


prm(source) = 1 
label(source) = 0.0 
curnod = source 
curtim 0.0 


C 
C Find the time of travel from the current node to each of the 

cS adjacent nodes which are not permanently Labeled. Find the minimum 
C distanced node, label it permanent. 

Cc 


do 1001 iter = 1, numnod - 1 


c print *, ‘iteration ', iter, ‘current node ', curnod, ‘time ', 
Cc & curtim 
c print *, 'minlbl ', minlbl 


minlbl = 99999 
C read (*,1) 2 
do 1002 node = 1, numnod 
if (prm(mode).eq.0) then 
if (curtim.ge.clabl(curnod, node)) then 
label (node) 


& = min(label(node), curtim + dtabl(curnod, node) ) 
else 
label (node) 
& = min(label(mode), ctabl(curnod,node) + dtabl(curnod,node) } 
endif 


if (label(node).lt.minlbl) then 
minlbl = label (node) 
minnod = node 


endif 
endif 
ic 
C 
1002 continue 
C 


curtim = label(minnod) 
prn(minnod) = 1 
curnod = minnod 
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do 9999 i = 1, numnod 
print *, i, label(i) 

continue 

read (*,1) gz 


QQqQaqg=maqaaAaAN 
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001 continue 
do 90 i = 1, numnod 
print *, 'node ', i, ‘prm(node) ', prm(i) 
0 continue ( 
return 
end 


C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C 
subroutine RTABLE(ctabl, dtabl, tail, head, c, d, atakd, 


& cdf, numtyp, numnod, numarc) 


parameter(maxnod = 40, maxarc = 60, maxtyp = 5) 


real ctabl(maxnod, maxnod), dtabl(maxnod, maxned), 
& c(maxarc, maxtyp), d(maxarc, O:maxtyp), cdf(maxarc, maxtyp), 


& rand 
integer numnod, numarc, numtyp, tail(maxuarc), head(maxarc), 


& atakd(maxarc), type 
common last 


Debugging print 
print. *, forvarce2-a. 
do 19 i = 1, numtyp 
write(*,2) i, -el(2, i), d(2, ye cdt(2,. 
9 continue 
format(’ destr type ', i3, ' e = *, £8.3,' d = ',f8.3, cdtie—uam 
& £8'.3:) 
read (*, 1) 2 


WOO OeCeaa a 


Initialize the table 


Le ® ay 


do 1001 i = 1, numnod 
do 1002 j = 1, numnod 
ctabl(i, j) = 99999 
dtabl(i, j) = 99999 
1002 continue 
1001 continue 


C 

C Assign the proper construction and time-of-travels in the tables 
C ctabl and dtabl 
C 


do 1005 i = 1, numarc 
if (atakd(i).eq.0) then 
C Assign the unharassed travel distance and no POTS: time 


dtabl(head(i), tail(i)) = d(i, 0) 
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ctabl(head(i), tail(i)) = 
dtabl(tail(i), head(i)) = d(i, Q) 
ctabl(tail(i), head(i)) = 
else 
p = rand({) 
C Debug print 
C 
C print *, ' p= ‘3 Pp 
C 
C Determine the type of destruction from the random probability 
C p, and assign the construction and travel times accordingly 
C 
do 1003 j = 1, numtyp 
if (p.it.cdf(i, j)) then 
type = 73 
goto 1004 
endif 
1003 continue 
1004 continue 
dtabl(head(i), tail(i)) = d(i, type) 
ctabl(head(i), tail(i)) = c(i, type) 
dtabl(tail(i), head(i)) = d(i, type) 
ctabl(tail(i), head(i)) = cli, type) 
endif 
1005 continue 
C 
C Debug print 
C 
C print *, 'Dtable and Ctable at the end of Rtable’ 
c do 1009 i = 1, numnod 
C write(*, 10) (dtabl(i, j), Jj = 1, numnod) 
COUY continue 
C write(*, *) 
C do 1008 i = 1, numnod 
c write(*#, 10) (ctabl(i, Jj), j = 1, mumnod) 
COQ08 continue 
G read (*,1) 2z 
return 
1 format(al ) 
10 Formac FOX fO.3, 2x), 7 FOCLE.3, 2x) 7 TOUiI8.3, 2x) / 
& 10(£8.3, 2x)) 
end 
C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C 
subroutine READNET(fname, head, tail, c, d, numtyp, numnod, 
& hnumarc, atakd, cdf, source, dest) 
parameter(maxnod = 40, maxarc = GO, maxtyp = 5) 


real c(maxarc, maxtyp), d(maxarc, O:maxtyp), cdf(maxarc, maxtyp), 
& prob(maxarc ) 
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integer numnod, numarc, numtyp, tail(maxarc), head(maxarc), 
& atakd(maxarc), arc, type, source, dest 

character*20 fname 

character*80 chrlin 


open(unit = 20, file = fname) 


Read the network data amd the attack plan 


OG 


read(20, 100) chrlin 

read(20, *) source, dest 
C print *, "s = °, Source, "t= deat 

read(20, *) numnod, numarc, numtyp 

read(20, 100) chrlin 
C print *, numnod, numarc, numtyp 
100 format(a80) 

do 1000 i = 1, numarc 

read(20, *) arc, tail(i), head(i), (c(i, type), type = 1, 


& numtyp), (d(i, type), type = 0, numtyp), atakd(i) 
1000 continue 
Cc 
c Read the probability of damage of each type and compute the cdf 
C 


read(20, 100) chrlin 
CG print *, ’numarc = °, numarc 
do 1001 i = 1, numarc 
read(20, *) arc, (prob(type), type = 1, numtyp) 
Cc print *, i, arc, (prob(type), type = 1, numtyp) 
edfi(i,*!) = sprob(!) 
do 1002 j = 2, numtyp 
edf(i, Jj) = edfC(l, J o— Py fapreb 


L002 continue 
1001 continue 
return 
end 
C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC CCG 
C 
subroutine outpt(arrive, numnod) 


C 
parameter(maxnod = 40, maxarc = 60, maxtyp = 5) 
C 
real arrive(maxnod) 
print * 
print x, Yess ees sst teste see Sess ees se ees sss sss sss sssssssssscsscs 
do 1000 i = 1, numnod 
print *, 'NODE ', i, ° time from source ', arrive(i) 
1000 continue 
return 
end 
C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
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Qa 


real function rand() 


last is the last integer generated by t 
parameter (L = 29, C = 217, xm = 1024.0) 
real p 
common last 


last = mod(last * 1 + c, xm) 
rand = real(last)/xm 


return 
end 
d 
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APPENDIX C. TABLES AND NETWORK DATA FILE (EXAMPLE 1) 


TABLE 4. RECOVERY TIME DATA (KOREA) 


fe [tone foo fo [ep To [eo [oe [ee 


road 
road 
bridge 
bridge 
mtn.road 
mtn.road 
road 
mtn.road 
mtn.road 
bridge 
mtn.road 
tunnel 
mtn.road 
road 
mtn.road 
bridge 
bridge 
road 
road 
bridge 
mtn.road 
bridge 
mtn.road 
bridge 
mtn.road 
road 
road 
mtn.road 
road 
road 
bridge 
road 
bridge 
road 
road 
road 
bridge 
road 
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TABLE 5. TRANSIT TIME DATA (KOREA) 


i) 13 1h 
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TABLE 6. NETWORK DATA FILE 


source, dest followed by numnod, numarc, numtyp 

1 26 

26 38 4 

arc, T, H, e(1),c(2),c3); e(4), d(0), d(1), 6 ° ’ 
1 1, 2; 0, O, 146.4, 244, 12.5,> Bi. 75, 1725, 
2 1, 4, O, O, 168.0, 280, 112.5, 223775,157.5, 
3 7a ts O, 0, 192.0, 320, 30.0, Wsa700, 42.0, 
4 2, 6; 0, O, 204.0, 340, 70.0, 7200, 98.0, 
5 Bess O, 0, 252.0, 7420; 40.0, 44.00, 56.0, 
6 a, 5, O, O, 237.6, 396, 32-5, Geo. to, 45.5, 
7 a, 12, O, OO, 184.8, 308, 45.0, 49.650, 63.0, 
8 Oni 6, 0, OO, 240.0, 400, 42.5, 36.76, 59.5, 
9 St afar 0, OO, 244.8, 408, 22.5, @eaeto, 315, 
10 (Br pctcae 0, 0, 216.0, 360, 42.5, "46.75, 59-5, 
11 ee 0, #0, 223-2. 312, 30.0, 33.00, 42.0, 
12 Un ln 0, “0, 732.0, 122¢. 52.5, @iv.m5, 73-25, 
V3). V2, is. 0, 0, 241.2, 8402, 42.5, 46.75, 59.5, 
14 9, 13, O, -0, 15122,.-252. 32.5, W5.05, 45.5, 
15 SS ae (oe O, O, 224.4,) 3a 35.0); 6.50, 4920, 
16 8 107 0, O, 219:6, 366, 55.0, 0.50, 77.0, 
LT. 27-08, O, 0, 189.6, 316, 60.0, G6.00, 84.0, 
18 13, 14, 0, O, 174.0, 290, 35.0, “38.50, #9.0, 
19 14, 15, O, -0,°163.2,.:272, 10.0, @ieCco,. 14.0, 
20 18, 19, O, O, 184.8, 308, 32.5, see7o, 45.5, 
21 14,°-19), 0, 0, 237.6, 395. 25.0, 2750, 35.0, 
22 5,516, O, OO, 192.0, 320. 50.05 65°00, 70.20, 
Z3. 10,716) 0, 0, 228.0,. 36¢, 52.5, 57:75, 73.5, 
24) 1857726, O, O, 194.4, 324, 119.55 lel.45,16723, 
25 IO y2)., 0, ©, 242.4, 404, 42.5, 46.75.. 59.5. 
26°) 1G-e-cOo. 0, 0,..150.0,--2507, 25.0, geet o0 35.0, 
27 VG? , ©, 0, 1523.6, 256. 37.5, 4025, 52.5, 
26) 10,1, 07" -0,°°252.05, 24207, 45.0, 49.50, 63.0, 
29 11, 17, O, O, 166.8, 278, 20.0, 42-00, 28:0, 
30°) 17,20, 0, 0, 163722. 272, 30.0, “35.00, 42.0, 
31) 2ie25 O, O, 189.6, 316, 32.5, WS .%o, 45.5, 
32 one lsreco 0, ©, 159-6. 266: 15.0, IG 750, 21.0, 
33 20,522, 0, 0, 198.0, 330, 62.5, “68.575, 87.5, 
34.) 22524. O, OO, 172.8, 280, 42.5, 46.75, 59.5, 
S078 225020, O, oO, 180.0, 3Ce, 12.5, “te.t5, 176, 
36 23, 24, 0O,. 0, 169.2,. 282, 30.0, 5.00, 4270, 
3 hme eae eo; O,. 0, 194.44, 324. 12.5, oe 75.,) 1720, 
38 25, 26, 0 O, 150.0, 250, 82.0, ‘90.20, 11476, 
are, poly p02); pert 4) 

1022 7580..2 0.550. 

2> (O22 uO e205 0 

Sr Oe Oe es, Ov Oye 

4. O03 27a0.2,. 0755.20.51 

5 “022,62, °0.5.,-0-1 

G6 0.20.2, °0.5,708! 

Tt 50.230. 2,2 025, 0 2! 

8°. 0.2,90.2,. 0,5, 001 

9. -0.:2,-0.2,-0-5. 0.4 
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d(4), atk (0 = 
13.75, 4.375), 

123475, 129.375. 
33.00, 34.5007 
77.00, 80.500, 
44.00, 46.000, 
35.75, -37. 35. 
49.50, S1.750, 
46.75, 48.875, 
24.75, 25.8752 
46.75, 48.875, 
33.00, 34.5002 
57.75, 60.3758 
46.75, 48.8752 
35.75, 87.3758 
38.50, 40.250, 
60.50, 63.250, 
66.00, 69.000, 
38.50, 40.250, 
T1006, 11.500, 
35.75, 37.3708 
27.50, 28. 750e 
55.00, 57.500, 
57.75, 60.375, 
131.45 137.4255 
46.75, 48.875, 
27.50, 28.7508 
41.25, 43.1255 
49.50, 51.7503 
22.00, 23.000, 
23.00, 34.5007 
35.75, 37.3750 
16.50, 17.2506; 
68.75, 71.875, 
46.75, 48.875, 
13.75, 14-3735 
33.00, 34.5005 
13.75, 14.375, 
90.20, 94.300, 


no 


oooo0o0coe°0c°oe°0c°009000o00o0 0 000O0OO0OOOFOOOOCOOKF OOF KF Ow 


oooooqoqooo°e°o°e°e°o 


> oe oe ©® 
039 09 09 8 DM OH HO ND OD DA DO DS 09 


O57. 1 
0.5, Ont 
OER) oO rag | 
On5, wet 
O.95;, 0.) 
0.5, 0.1 
Of 5 5) Oar 
O75, O21 
0.5; Oe 
0.6, 0.1 
GO. 5, 0.1 
0.5.0.) 
O45, O.1 
On oy Cat 
O.5., Ore 
0.5, 0.1 
0..5,°0.1 
O70, O.1 
O...5:- 2 Oar 
O25, O01 
Oc5, Ol 
G.5, O.1 
O75, Oca 
O25, 0.1 
Ovo, O12 
0.5, 0.1 
O54) 0.1 
025. O74 
O25, Oct 
* c(4) 
* d(0) 
* d(0) 
* d(0) 
* d(0) 


TABLE 6 (cont.) 


51 


TABLE 7. TRANSIT TIME COMPARISON BETWEEN TWO MODELS 
(Scenario 1) 


Determin. Stochastic 
Arc Set 
Time C195 %) | C.1(90 %) 

42 b6nlen 346.9000 368.6438 (366.6758, (366.9921, 
370.6118) 370.2955) 

24,25, 26,28 356.7175 360.9398 (358.2446, (358.6778, 
363.6351 ) 363.2019) 

23 Ole 377.0000 371.3937 (369.8125, (370.0666, 
372.9750) 372.7209) 

PS ITE Bhs: 418.0300 401.7033 (399.0912, (399.5110, 
404.3153) 403.8955) 

BESTE AU 414.5000 401.2283 (399.0602, (399.4086, 
403.3964) 403.0479) 

2, 4,10,38 427.0000 391.5468 (387.8984, (388.4847, 
395.1952) 394.6089) 

Z, 410316 467.5625 404.2913 (397.2708, (398.3990, 
411.3118) 410.1835) 

2 old 467.5625 416.2100 (407.2977, (408.7300, 
425.1223) 423.6899) 

ae. BAG 414.5000 393.2225 (389.7153, (390.2789, 
396.7297) 396.1661) 

TAG 355.1250 379.3322 (375.0317, (375.7228, 
383.6328) 382.9416) 
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TABLE 8. TRANSIT TIME COMPARISON BETWEEN TWO MODELS 
(Scenario 2) 























Arc Set 
Time C195 %) | C.1(90 %) 
4,12,16,17 | 397.0000 | 378.9950 | (376.0747, | (376.5440, 
381.9153) 381.4460) 
24,25,26,28 | 398.4175 | 399.5087 | (394.6451, ‘| (395.4267, 
404.3723) | 403.5907) 
2, 3, 6,12 | 377.0000 | 371.3937 | (369.8125, | (370.0666, 
372.9750) 372.7209 
2, 424,38 | 418.0300 | 419.2902 | (414.5211, ‘| (415.2875, 
424.0593) 493.2928) 
2,4, 7,10 | 414.5000 | 401.2563 | (399.0849, | (399.4338, 
403.4276) 403.0787) 
2, 410,38 | 427.0000 | 391.5508 | (387.9015, _ | (388.4880, 
395.2000) 394.6136) 
2, 4,10,16 | 579.0625 | 427.0912 | (412.2838, | (414.6636, 
441.8987) 439.5189) 
2,4, 5,10 | 579.5625 | 454.6700 | (435.9646, | (438.9708, 
473.3754) 470.3692) 
2, 4,812 | 414.5000 | 393.2225 | (389.7153, | (390.2789, 
396.7297) 396.1661) 
7,11,12,16 | 478.3250 | 405.5908 | (393.2462, | (395.2302, 
417.9355) 415.9515) 
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APPENDIX D. TABLES AND NETWORK DATA FILE (EXAMPLE 2) 


TABLE 9. RECOVERY TIME DATA (CALIFORNIA) 


fe [tomo fo [een Tap Teo [ome [= 


Co 1MHo LR WNrF 


— — rt TD DDI DOD RS —— — 
a 


— 
prod 
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NOVMWNVUNNONY NH NHN NH H 


S290000 CCC GCa0 


TABLE 11. 


NETWORK DATA FILE 


dest followed by numnod, 


arc,p(1),p(2),+-; 


H, c(1), 
2; 0, 
3 0, 
35 0, 
4, 0, 
O's 0, 
6, 0, 
6, 0, 
7, 0, 
1; 0, 
8, 0, 
9, 0, 
sie 0, 
10, 0, 
sig Gee 0, 
ie 0 
: <2 
O22 420%.-5, 
Q.2,.0.5, 
O22, 0.5, 
On 2 080 ; 
0.2, 20% 5:, 
Orc, 0.5, 
Ome, 0.5, 
Ole, O15: 
Ome, 0.5, 
0) A oe 
Orc. 0.0, 
One) 5 
Oey Os), 
Owe. 0.5; 
0 * (4 
1.1 * d(0 
1.4 * da(0O 
1 * d(0 
1 * d( 


— 


Koo Ro Ro Roo Ro Ro Rohe Nolo t— i 
ee ee or onan 


© weer wee eee we 
~— 


c(4), 


190, 
252, 
204, 
260, 
260, 
190, 
246, 
258, 
252, 
306, 
326, 
266, 
270, 
254, 
256, 


numarc, 


d(0), 


76, 


dO 


numtyp 
GG), ai. (4)y atk 
3.6, 106.4, 63.6, 
39.6, 50.4, 39.6, 
61.6, 78.4, 61.6, 
8.8, ec 8.8, 
39.6, 50.4, 39.6, 
2 ao GT. 2 ; O22 ac; 
2e.0, 2520), 2c. 
-2, 44.8, BS Ac 
39.6, 50.4, 39.6, 
66.0, 84.0, 66.0, 
SO 56.0% 44.0, 
6, 90.4, 0. 39.6, 
16, eecee, iG, 
SOOO. sc, LOZec, 
29, lo2aoeeel lS .9, 


~=—— 


COCO COOFFFH DQG O0O00 


TABLE 12. TRANSIT TIME COMPARISON BETWEEN TWO MODELS 
(Scenario 3) 


Stochastic 


(337.0615, 
351.0530) 


ee! 367.0850 346.1183 (337.0435, (338.5019, 
355.1932) 353.7347) 

357.6800 346.9688 (338.6827, (340.0144, 
355.2550) 353.9233) 


2,3, 7 309.0000 279.2960 (276.6637 (277.0867, 
281.9283) 281.5053) 

1,2, 8 357.0200 338.7611 (329.9167, (331.3381, 
347.6056) 346.1841) 

ob 358.3400 346.1733 (337.6795, (339.0446, 
354.6671) 353.3020) 

fo 273.4800 282.9700 (279.5856, (280.1295, 
286.3544) 285.8105) 
2,5, 8 293.0000 291.0300 (290.1343, (290.2783, 
291.9257) 291.7817) 
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TABLE 13. TRANSIT TIME COMPARISON BETWEEN TWO MODELS 
(Scenario 4) 


Arc Set 

8, 9, 10 376.5400 323.8200 (Bible 1291, (313: Wereifee 
336. 5109) 334. 4713) 

Pe 2, 452.5400 397.9328 (379.6254, (382.5677, 
416. 2402) 413. 2979) 

iz, 4 453.8600 401.3493 (384. 2001, (386. 9562, 
418. 4986) Ao 7424) 

lee, 14 467.8850 409.0539 (391.7082, (394.4959, 
426. 3995) 423. 6118) 

457.5400 403.1413 (386.1323, (388.8659, 
420. 1503) ats 4167) 

Be O, 293.0000 279.2960 (276.6637, (277.0867, 
281.9283) 281.5053) 

m2 8 457.5400 395.0933 (377.4376, (380.2751, 
412.7490) 409.9115) 

lhe 2.11 459.1400 404.1253 (386.7449, (389.5382, 
421.5057) 418.7124) 

a8. 9 309.0000 282.9700 (279.5856, (280.1295, 
286.3544) 285.8105) 

By o, 8 293.0000 291.0300 (290.1343, (290.2783, 
291.9257) 291.7817) 
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